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DECIGO Pathfinder (DPF) has an ability to detect gravitational waves from galactic 
intermediate-mass black hole binaries. If the signal is detected, it would be possible to de- 
termine parameters of the binary components. Furthermore, by using future space-borne 
gravitational wave interferometers, it would be possible to test alternative theories of 
gravity in the strong field regime. In this review article, we first explain how the de- 
tectors like DPF and DECIGO/BBO work and discuss the expected event rates. Then, 
we review how the observed gravitational waveforms from precessing compact binaries 
with slightly eccentric orbits can be calculated both in general relativity and in alter- 
native theories of gravity. For the latter, we focus on Brans-Dicke and massive gravity 
theories. After reviewing these theories, we show the results of the parameter estimation 
with DPF using the Fisher analysis. We also discuss a possible joint search of DPF and 
ground-based interferometers. Then, we show the results of testing alternative theories of 
gravity using future space-borne interferometers. DECIGO/BBO would be able to place 
4-5 orders of magnitude stronger constraint on Brans-Dicke theory than the solar system 
experiment. This is still 1-2 orders of magnitude stronger than the future solar system 
mission such as ASTROD I. On the other hand, LISA should be able to put 4 orders of 
magnitude more stringent constraint on the mass of the graviton than the current solar 
system bound. DPF may be able to place comparable constraint on the massive gravity 
theories as the solar system bound. We also discuss the prospects of using eLISA and 
ASTROD-GW in testing alternative theories of gravity. The bounds using eLISA are 
similar to the LISA ones, but ASTROD-GW performs the best in constraining massive 
gravity theories among all the gravitational wave detectors considered in this article. 
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1. Introduction 

1.1. Gravitational Waves in General 

In general relativity (GR), gravitational fluctuations propagate as "ripples" in the 
spacetime at the speed of light, known as gravitational waves (GWs). Since gravity 
is weak compared to other fundamental forces, the interaction of GWs with other 
matters is very faint. This means that GWs can escape from highly dense or optically 
thick regions that cannot be probed by electromagnetic (EM) waves (e.g. very early 
universe before the last scattering surface). Therefore, GWs have potentials to open 
up a novel astronomy and cosmology. 

The existence of GWs has been proved indirectly by measuring the orbital decay 
rate of binary pulsars, caused by the energy loss through gravitational radiation. 
This can be calculated analytically within the framework of GR, and this matches 
with the measurements beautifully^'^. However, this indirect GW measurement only 
detects the radiated energy via GWs and not the perturbation of the spacetime 
itself. Therefore, the direct detection of GWs has been awaited for a long time to 
detect this latter effect. Currently, several ground-based detectors that are aiming 
for the direct detection are shifting from initial phases (flrst generation) to advanced 
phases (second-generation) such as adv. LIGO^-'', adv. VIRGO^'^, GEO-HF'''® and 
KAGRA (formerly called LCGT)^'^°. These are GW interferometers that consist of 
two arms. When GWs pass through this detector, the length of each arm changes, 
modifying the interferometric pattern. KAGRA is novel and sometimes called 2.5th- 
generation detector since it will be cryogenically-cooled and is buried underground 
to reduce the thermal and seismic noises. These detectors have their best sensitivities 
at around 100-lOOOHz. Basically, the high-frequency parts of these detectors are 
limited by the shot noises while the low-frequency parts are limited by the radiation 
pressure and the seismic noises. 

Sources and sciences from GW interferometers have been summarized in 
Refs."'^"'^' Promising sources for these second-generation detectors are signals 
from compact binaries (see Sees. 3.4 and 3.5 of Rcf.^^ for a theoretical overview 
and Refs."'^^"^^ for observational results) whose event rates are estimated as 
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C'(10)/year-'^''. (See Ref.^^ for the review on the electro-magnetic counterparts of 
GW sources.) Binary signals are divided into three phases, inspiral, merger and 
ringdown. During the inspiral phase, the binary separation gradually shrinks due 
to the energy loss via GW radiation. The inspiral waveform has been studied well 
under the post-Newtonian (PN) formalism Merger phase is highly non-linear 
but there has been a great progress thanks to numerical relativity^^^^** . Ringdown 
phases are also well-studied using black hole (BH) perturbation method^^. GWs 
from a merging compact binary is often called chirp signals since their frequencies 
get higher and higher as they head towards merger. Since the inspiral waveform 
modeling has been performed well using PN approach^^, it can be used to estimate 
binary parameters'^ and to probe cosmology as standard sirens''^. Also, neutron- 
star (NS) binaries can be used to constrain equation of state (EoS) of NSs'®"^''. 
Other interesting sources are burst signals (see Sec. 3.2 of Ref.^' for a theoretical 
review and Refs.'^^^^^ for the current observational results) including supernovae^', 
magnctars and cosmic strings''^' "^^j and continuous signals (see Sec. 3.3 of Ref.^' 
for a theoretical overview and Refs.^^^^^ for observational results) from e.g. (newly- 
born) neutron stars. There are also stochastic GW backgrounds (see Sec. 3.6 of 
Ref.^' for a theoretical review and Ref.*® for observational results) as potential 
candidates for GW sources such as the ones from cosmic-strings''^, pre-big-bang 
models^" and preheatings^'^. GW backgrounds associated with inflation may be de- 
tected by ground-based detectors if the spectrum has a sufficient blue tilt. There is 
also a proposed third-generation project called Einstein Telescope (ET)^'-^^. The 
sensitivity is increased roughly by one order of magnitude compared to the second- 
generation ones and ET has interesting sciences which have been summarized in 
Ref.54. 

1.2. Space-Borne Interferometers 

Ground-based interferometers have difficulty in detecting GW signals lower than 
1-lOHz due to the seismic noise. To overcome this problem, space-borne interfer- 
ometers have been proposed. Among them, the Laser Interferometer Space Antenna 
(LISA) has been proposed by ESA and NASA with its optimal sensitivity around 
GW frequency of ImHz^^'^^. Classic LISA consists of three spacecrafts from which 
laser is emitted to one another. This triangular-shaped interferometer has arm- 
lengths of 5 X lO^kni. Due to this extremely long armlength, LISA is a transponder- 
type (rather than a Michelson or a Fabry-Perot-type) interferometer. (I.e. Instead 
of reflecting the incoming laser directly, each spacecraft once receives it and emits 
the laser with corresponding phase.) The triangular cluster follows the same orbit 
as the Earth, keeping its position 20° behind the latter. The expected targets are 
supcrmassivc black hole (SMBH) binaries and white dwarf (WD) binaries. The for- 
mer detection may help in clarifying the mechanism of how SMBHs are formed. A 
system consisting of a small compact object orbiting a SMBH is called an extreme 
mass ratio inspiral (EMRI)^*", and GWs from this system encodes the information 
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Fig. 1. The configuration of DECIGO and BBO**. There are eight efl'ective interferometers in 
total. 

of the structure of the strong-field spacetime^®"^^ . Unfortunately, GWs from WD 
binaries may mask*"^' ''^ other signals including primordial GW background''^' As 
a prototype mission, LISA path finder (LPF) will be launched before LISA®^. 

Although, the classic LISA was planned and developed by NASA and ESA, the 
collaboration between these 2 agencies terminated in 2011. The new version of LISA 
(called eLISA) led by ESA alone is now under consideration. The arm-length would 
be shortened by 5 times compared to the classic one and the number of the arms 
would be reduced to 2. The satellites are slowly drifting away from Earth to save 
propellant. The expected operation period is 2-5 yrs. See Refs.®^"^'' for expected 
sources and sciences with eLISA. 

Another space-borne interferometer, the Deci-Hertz Interferometer Gravita- 
tional Wave Observatory (DECIGO) has been proposed by Seto, Kawamura and 
Nakamura®^. It consists of four clusters, each having a triangular interferometer. 
The major differences between DECIGO and LISA are that the armlengths are 
shortened to lO'^km and the former is a Fabry-Perot-type interferometer whereas 
the latter is a transponder-type as previously mentioned. Two out of four clusters 
are situated on the same site so that the correlation analysis can be performed to 
detect primordial GW background, while the other two are placed far apart so as 
to increase the angular resolutions of the sources^""^^ (see Fig. 1). It is most sensi- 
tive at around O.l-lHz. Similar interferometer, the Big Bang Observatory (BBO) 
has been suggested as a follow-on mission to LISA^'^'^^. Since it is a transponder- 
type interferometer, it can be considered as a mini-LISA with its configuration 
similar to the one of DECIGO. (The noise curves of adv. LIGO, ET, LISA and DE- 
CIGO/BBO are shown in Fig. 2 ^.) Since the WD/WD binary signals has a high 
cutoff frequency at / ~ 0.2IIz^^, the main target of DECIGO/BBO is the primordial 



^DECIGO has been proposed with a few times less sensitivity than BBO, but this is not a fixed 
design and here, we have assumed that it has the same sensitivity as BBO. 
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Fig. 2. The square root of the noise spectral densities Sn{f) of various detectors. (The noise 
spectral density roughly corresponds to the lower limit of the GW amplitude that each interfer- 
ometer can detect.) The (blue) thick dashed and the (red) thick solid curves represent the LISA 
and DECIGO/BBO ones, while the (green) thin solid and the (purple) thin dashed curves show 
the ones for ET and adv. LIGO, respectively. The (orange) dotted-dashed and the (black) solid 
curves correspond to the cLISA and ASTROD-GW ones, respectively. 



GW background (PGWB)^^'^^'^^. NS/NS foreground signals might mask PGWB 
instead of WD/WD signals, but it is likely that BBO would have enough sensitivity 
to subtract sufficient amount^^'^®. See Refs.^^'^''''^^"^-'^* for other cosmological and 
astrophysical scientific potentials of DECIGO and BBO. 

As for the prototype mission, DECIGO Pathfinder (DPF)^!^ is hoped to 

be launched in 2016-2017. The main goals are to test the key technologies for the 
space mission and to carry out observations of GWs and Earth gravity. The first 
space-borne GW detector, SWIM, has already been launched in 2009 n^, ii7^ 
is a torsion-bar type space antenna. The ground-based torsion-bar antennae have 
been used to place upper limits on the energy density of stochastic GWs at O.l-lHz 
118, 119^ The improved version of the torsion-bar antenna (TOBA) has been proposed 
by Ando et al}^"^ . 

DPP has a sensitivity of y/SJJ) - 2 x lO^^^ Hz"^/^ at / ~ 1-lOOHz. The 
main GW sources for DPF are the intermediate mass BH (IMBH) binaries. There 
are evidences that stellar-mass BHs and SMBHs exist, but there is no unambiguous 
detection of individual IMBHs. IMBH detection may reveal the formation mecha- 
nism of SMBHs. DPF has enough sensitivity to detect GWs from galactic IMBH 
binaries that might exist at the centers of globular clusters and massive young clus- 
ters (see Refs.^^°'-^^^ for reviews on IMBHs and Ref.^^^ for the possibility of IMBH 
binary formation). For the observation of the Earth gravity, DPF has an ability to 
perform complementary operation compared to other missions such as CHAMP^^"^, 
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GRACE124 and GOCEi^^. 

ASTROD-GW (Astrodynamical Space Test of Relativity using Optical Devices 
optimized for GW detection) is another space-based GW interferometer mission. 
It consists of 3 satellites that are placed at near Lagrange points L3, L4 and L5, 
respectively, forming a large triangluar interferometer similar to LISA. Since its 
armlength is 52 times larger than that of LISA, it is more sensitive than LISA at 
lower frequency. 

1.3. Testing Alternative Theories of Gravity using GWs 

Among all the sciences that have been mentioned above, there is yet another very 
interesting and fundamentally important science that can be best probed by using 
GWs. It is the tests on alternative theories of gravity^^'^^^^^ . In order to solve prob- 
lems like dark energy-'^^'', dark matter and inflation^'^-'^' ^■^^ within the context of GR, 
usually we need to introduce unknown matters or fields, but there are possibilities 
that these problems can be explained naturally by modifying gravitational theories. 
Also, if the classical gravitational theory is to be realized at the low energy limit of 
more fundamental theory like superstring theory ^'^^^ the classical one does not 
necessarily reduce to GR due to the additional fields (e.g. dilatons) that couple to 
gravity. 

The tests on alternative theories of gravity have been performed with great 
accuracies under the solar system experiments and binary pulsar observations^' 
and no deviation from GR has been reported so far. The former can probe only the 
weak field limit of the theory, while the latter can probe the strong field gravity^'^^ 
such as the effacement principle and the scalarization^'^^"^^® in scalar-tensor theory. 
However, since the typical velocity of the binary component is v/c ~ 10""^, the 
system is not so dynamical. We are here interested in testing GR in the strong and 
dynamical field regime. The best way to perform these tests in the near future is 
to use GWs (especially GWs from compact binaries) since they directly contain 
gravitational information in the strong-field regions. Also, GWs lead to remarkably 
high precision GR tests in the following sense. Usually, GWs are so weak that they 
are buried under the detector noises. In order to dig them out from the noises, 
we need to perform the matched filtering analysis^^'^^^ . This is nothing but taking 
correlations between the observed signals and the template waveforms so that if 
the former match with the latter, we can extract the signals out of the noises. 
Therefore, parameter estimations using this matched filtering technique are sensitive 
to deviations in GW phase since the correlations are remarkably reduced even if the 
phases between GWs and the templates are different only slightly. Now, it is often 
the case that alternative theories of gravity modify the phases of GWs from compact 
binaries (see Fig. 3). Let us say that we observe a signal at a frequency / — O.lHz 
for observation period tobs = lyr- Then, the number of GW cycles (the number of 
phases) that wc see is roughly / x iobs ~ 3 x 10^. This shows that we would be able 
to detect the deviations from GR if they modify the phase at least by (3 x 10^)^^ ~ 
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time 



Fig. 3. A schematic picture of gravitational waveforms from a compact binary in GR (red thick) 
and in alternative theories of gravity (blue thin). 



3 X 10~^! There are many works calculating how accurately we can probe alternative 
theories of gravity using GWs. (See e.g. Ref.^^° and references therein, and see 
Refs.^'"'"^^^ for model- independent tests. See Refs.^^'^"^^^ for possible tests on GR 
using ringdown signals.) 

Two of the most important characteristics of GR are (i) it has only two tensor 
gravitational degrees of freedom and (ii) the graviton is massless (in other words, 
GWs propagate at the speed of light). In this review, we focus on simple extensions 
of each of the above points by (i) introducing an additional scalar field and (ii) 
giving a finite mass to graviton. 

As for the first point, we consider Brans-Dicke (BD) theory^^® which is the 
simplest representative of the scalar-tensor theory This theory is parameterized 
by the so-called Brans-Dicke parameter wbd which gives the (inverse) coupling 
between the additional (massless) scalar field and the matter fields. If we take the 
decoupling limit of ojbd ^ oo, this theory reduces to GR. The current strongest 
constraint has been obtained from the solar system experiment by Shapiro time 
delay measurements using the Cassini satellite^^° as wbd > 4 x 10^. There is a 
future space mission called ASTROD which is expected to be launched in 2021^^^. 
It would perform very precise tests of GR in the solar system and the constraint 
on ojbd would be 3 orders of magnitude stronger than the Cassini bound. One 
of the remarkable points about considering gravitational radiation in this theory 
is that due to the additional scalar field, there exists scalar dipole radiation^^^"^^''^! 
Because of this additional energy release, the binary evolution is modified compared 
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to GR; which further makes a change in GW phase. The orbital decay of the binary 
pulsar PSR J1738-0333 has put a slightly weaker constraint than the solar system 
one^^^. There have been several works calculating the possible bounds on wbd using 
future GW interferometers^^^"^™. (See Refs.^^^' ^''^ for the constraint on massive BD 
theory, Ref.^^^^ for the one on generic scalar-tensor theories and Ref.^^* for the bound 
on generic dipolar gravitational radiation.) Unfortunately, these results show that 
the bounds from future GW observations using adv.LIGO and LISA would not be 
able to exceed the one from the solar system experiment. In this review, we will 
show that DECIGO/BBO would be able to place 4 orders of magnitude stronger 
constraint than the Cassini bound"'^''^. 

For the second point, there exist many theories in which graviton acquires a 
finite mass nig i'^6-i82^ ^g^g ^ g_ j^^^g I83,i84 recent reviews.) These theories are 
called massive gravity or massive graviton (MG) theories. From a pioneering work 
by Fierz and Pauli^^^ in 1939, most of massive gravity theories have suffered from 
problems like Boulware-Deser ghost modes appearing in the curved background^®^ 
or violating the Lorentz invariance^^*' Only recently a self-consistent massive 
gravity is proposed^^^' Independent of these specific massive gravity theories, 
the gravitational potential is modified to the one of Yukawa type. This means that 
the effective gravitational constant depends on the distance from the source, which 
modifies the Kepler's third law. This has been tested in the solar system and the 
constraint has been put on the graviton Compton wavelength as \g > 2.8 x 10^''cm 
As for GWs, the propagation speed is modified from the speed of light which 
again makes a change in GW phase. GWs can also perform a model-independent 
tests on the graviton mass. Compared to the solar system experiment, ground-based 
detectors can only put comparable (or slightly stronger) constraints^*''""'^*^ whereas 
LISA bounds should be stronger by four orders of magnitude^'*^' ^^''""'^^^! 
This is because the deviation in the propagation speed of GWs is larger for lower 
frequency GWs. DECIGO/BBO would be able to place slightly weaker constraint 
than LISA^^^ and DPF may be able to place comparable constraint to the solar 
system bound 



1.4. Organization 

In this review, we will first explain how accurately DPF would be able to determine 
binary parameters by detecting GWs from galactic IMBH binaries. Then, we will 
describe the ability of testing GR using space-borne GW interferometers, especially 
DECIGO/BBO and DPF. This article is organized as follows. In Sees. 2 and 3, we ex- 
plain the basic concepts, noise curves and event rates of DPF and DECIGO/BBO, 
respectively. In Sec. 4, we derive GWs from compact binaries, starting from the 
quadrupole, leading contribution for the binary with a circular orbit, and extending 
it to include higher PN terms and a slightly eccentric orbit. Then, in Sec. 5, we ex- 
plain how to construct the observed gravitational waveforms for the spin-aligned and 
precessing binaries. We also introduce the inspiral-merger-ringdown hybrid wave- 
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form. In Sees. 6 and 7, we eonsider BD and MG theories, respeetivcly. We explain 
how the GW phase is modified from GR, and review eurrent constraints. In See. 8, 
we show the results of the binary parameter estimation using DPF, and in Sec. 9, we 
show the proposed constraints on BD and MG theories (and also on other theories) 
using future GW interferometers. Finally, we conclude in Sec. 10. 

2. DECIGO Pathfinder 

2.1. Basic Designs 

]3pp9i,ii5 jg prototype mission of DECIGO to test the advanced technologies 
of GW space mission such as (i) a precise position measuring system with Fabry- 
Perot (FP) cavity, (ii) a highly stabilized laser source and (iii) the drag-free control 
system which shields external forces due to solar radiation and residual atmosphere. 
The weight of the satellite is about 350kg and it will be orbiting the earth at the 
Keplerian velocity with an altitude of 500km. The expected observation period is 1 
year. It consists of 1-arm interferometer with armlength 30cm and the laser power of 
lOOmW. FP cavity has not been tested in space up to now and DPF is expected to 
have better sensitivity than LPF which uses a Mach-Zender interferometer^^. DPF 
provides new possibilities for a precise measurement of position and high-stabilized 
laser in space. 

Each mirror is placed inside a module called housing. The relative positions of 
these mirrors and the frame will be measured by the local sensor and will be fed back 
to the satellite position using thrusters. LPF will operate at the Lagrange 1 (LI) 
point where the gravitational environment is stable, while DPF will demonstrate it 
in an earth orbit. This will open a new window for future space missions. 

DPF GW observations aim at the frequency of O.l-lHz which is important be- 
cause the ground-based interferometers and LPF are not sensitive in this frequency 
range. Also DPF data are expected to be more complicated than the one from 
ground-based interferometers due to the satellite orbital motion and the effects of 
the earth, and hence the development of data analysis technique for DPF has a 
significant meaning. 

2.2. Noise Spectrum and Observable Range of DPF 

The detected signal s{t) — h{t) + n{t) can be expressed as the combination of the 
GW signal h{t) and the noise n{t). If the noise is stationary, the noise spectral 
density Sn{f) can be defined as 

{n*{mf))^5{f-n\snU). (1) 

where n{f) represents the Fourier component of the noise and the angle brackets 
denote the expectation value. The root noise spectral density Sn{f) of DPF is 
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Fig. 4. The root noise spectral density \/ Sn (/) of DPF is shown as the (red) thick solid curve, 
together with the one for adv. DPF (high frequency part is limited by the shot noise rather than the 
laser frequency noise) as the (green) thin solid curve. We also show the adv. LIGO low frequency 
cutoff lOHz as the (magenta) thick dotted line and the one for KAGRA down to IHz as (blue) 
thin dotted curve. The sky-averaged amplitude of a IO^Mq equal-mass binary in u) Centauri is 
depicted as (light blue) thick dotted-dash curve and the one of a 10'' Mq equal-mass binary in 
NGC 6752 is drawn as (black) thin dotted-dashed curve. Here, we assumed that the observations 
are carried out from / = 0.03Hz to mergers. We also assumed that these BHs have dimensionless 
effective spin angular momentum x = 0.2. (This figure is taken from Rcf. ^®*.) 



Sn{f)^^-^^^^(4r] + 4.0 X 10-3" Hz-\ (2) 



given as 



(27r)4 VlHz^ 

where the first term corresponds to the acceleration noises while the second term 
represents the laser frequency noise. It is shown in Fig. 4 as the (red) thick solid 
curve. We set the lower frequency cutoff at / = O.OSHz due to the Earth gravity. 
It may be possible to improve the sensitivity in which the higher frequency part is 
now limited by the shot noise. We call this the "adv. DPF" whose noise spectral 
density is given by^^^ 

SnU)^-'^^ = '-^^^ (l£) " + 9.0 X 10-- Hz-^ (3) 

It is shown in Fig. 4 as the (green) thin solid curve. We also show the lower cut- 
off frequencies of the adv. LIGO (magenta thick dotted) and KAGRA (blue thin 
dotted) ^. 

In Fig. 5, we show the (sky-averaged) observable range of DPF and adv. DPF 
as the (red) thick solid curve and the (green) thin solid curve, respectively. We also 



''It is likely that adv. LIGO and adv. VIRGO also have non-zero sensitivity at / < lOHz, but 
KAGRA seems to have better sensitivity at this frequency range. 
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Fig. 5. This figure sliows the observable range for each detector with SNR threshold pthr = 5 
The (red) thiclc solid curve represents the pattern-averaged one for DPF while the (green) 
thin solid curve shows the one with adv. DPF. Other 2 curves correspond to the same detectors 
as in Fig. 2 "*" shows the masses and distances for assumed equal-mass IMBH binaries in 
the galactic globular clusters shown in Table 1, while "x" shows the ones displayed in Table 2 of 
Ref.^^*, whose masses have been determined by using M — a relation. 



Table 1. The distances and possible IMBH masses 
at the centers of galactic globular clusters 



NGC 


distance 


(total) BH mass 


No. 


(kpc) 


(Mq) 


5139 (w Cen.) 


4.8 


(3.0-4.75) X 10* 1^' 






<1.2xl0* 198'"9 






(1.3-2.3) xlO* 200 


6388 


10.0 


5.7x10^ 202 


6715 (M54) 


26.8 201 


9.4x10^ 203 


6752 


4.0 201 


2.0x10^ 204 


7078 (M15) 


10.3 201 


3.2x10^ 205 



show the ones of adv. LIGO and KAGRA with the same curves as in Fig. 4. The 
(black) dashed horizontal line at Z?^ = 8.5kpc corresponds to the galactic center. 
We also plot the possible GW sources at the centers of the globular clusters shown 
in Table 1, assuming that they consist of equal- mass IMBH binaries. 

2.3. Beam Pattern Functions and the Effect of Detector Motion 

The beam pattern functions of 1-armcd interferometer are given by^^^ 

-F"^ = ^ sin^ 9 cos 2^ , (4) 
F"" {9, ^) = ^ sin^ 9 sin 2-4} . (5) 
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Fig. 6. Wo introduce the detector coordinates {x, y, z} with its origin placed at the center of the 
DPF arm. The x-axis coincides with the detector arm while z-axis is orthogonal to this x-axis and 
exists on N-x plane, with N denoting a unit vector pointing towards the source. The angle of 
N measured from x is denoted as 9 while the polarization angle ip is defined as the angle from 
p = ^^^"^ to n which is a unit vector made by projecting z onto the plane perpendicular to N. 
fi is orthogonal to N and lies on x-z plane. (This figure is taken from Ref.^®*.) 



Here, 9 is the angle between the arm and the incoming GW, and i}} represents 
the polarization angle (see Fig. 6). The sky-averaged values of them are {F^"^) ~ 
^^x2^ = 2/15. In Fig. 4, we also show the sky-averaged GW amplitudes of a 
lO*Af0 equal-mass IMBH binary in uj Centauri (light blue thick dotted-dashed) 
and a IQ^Mq equal-mass one in NGC 6752 (black thin dotted-dashed). We assume 
the dimensionless spin parameter x of these BHs to be x = 0-2 whose value is 
predicted to be the most probable^"^ (see also Ref.^°^). We use the phenomenological 
inspiral-merger-ringdown hybrid waveforms^°^'^°^ (see Sec. 5.3) to estimate these 
amplitudes. One can see that GW frequency of IMBH binary in a; Centauri is too 
low for the ground-based interferometers, and hence its GW signals become unique 
sources for DPF. On the other hand, the one in NGC 6752 can be detected with 
both DPF and the ground-based ones. This suggests that it may be possible to 
perform joint searches between these detectors^^^, which we explain in Sec. 8.1. 

For later use, we introduce the celestial coordinates {x, y, z\ shown in Fig. 7. 
ai-axis points the vernal equinox, while z-axis points the north celestial pole and is 
orthogonal to the celestial plane. The source direction and the polarization angle 
(0, '0) can be re-expressed in terms of the source direction (^Si 0s) a-nd the direction 
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ecliptic 



vernal 
equinox 



Fig. 7. Wc introduce the celestial coordinates {x, y, z} with S-axis pointing the vernal equinox 
and the S-axis is orthogonal to the celestial plane and is pointing the celestial north pole. DPF 
orbits the earth with its plane always facing towards the sun and its arm pointing the earth. It 
has an inclination angle «d = 97.3° against the celestial plane. The magnified picture of the DPF 
orbital plane is shown at the top left panel. The angular velocity of the earth lu^ is measured from 
the X-axis while the one of the DPF cjd is measured from the axis which is obtained by projecting 
the 2-axis onto the DPF orbital plane. (This figure is taken from Ref.^^*.) 



of the 
cos 



as 



194 



(6) 



orbital angular momentum 0l 
= cos ipB (t) sin cos 

- [cos{i^E(i) " 0s} cos(y9D(i) cosiD + sin{(pE(i) - (^s } sin (^d (<)] sin 0s 
cosV' = cos^s{— cos{ipE{t) — 0l} sin ip£,{t) + sm{(pE{t) — 4>l} cos(^d(0 cos^d} sin^L 
+{cos{i^E(i) - 4>s}sm(p-D{t) COS^L 

— COS (pY){t){sm{(pEit) — 0s} COS cosiD + sin{0L — 0s} sin^L sin id}} sin^s 
(1 — cos'^ 9){1 — (sin^s sin^L cos(0l — 0s) + cos^s 



jcos^l)'} 



-1/2 



(7) 



2.4. Event Rate 

In this subsection, wc estimate the event rate for (I) Equal-mass BH binaries and 
(II) IMRIs using DPF. 



2.4.1. Equal-mass BH Binaries 

DPF can detect GW signals from galactic IMBH binaries. Two of the candidates 
where IMBH binaries might exist in our Galaxy arc at the centers of globular clusters 
(GGs) and galactic massive young clusters (GMYGs). 
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For the former case, we select 21 globular clusters out of 150 globular clusters 
that have been found. We determine the total mass of the possible IMBH binary 
by applying the M — a formula obtained by Tremaine et al.'^^^ as 

^^ = 1-^^^ (200^)'^- 

We assume the mass ratio of IMBH binary to be 1. The results are listed in Table 2 
of Ref.^^*. We plot IMBH binary of each globular cluster in Fig. 5 and count how 
many of them lie within the observable range of DPF and Adv. DPF. One sees that 
2 out of 26 globular clusters (5 shown in Table 1 + 21 shown in Table 2 of Ref.^^^) 
might contain IMBH binaries detectable with DPF on average. Then, the detection 
rate of the IMBH binaries in globular clusters Ngc is given by 

Here, following Ref.^^^, we divide the number of globular clusters by the age of the 
universe. This is because only one IMBH binary can be formed over its lifetime in 
each cluster. When we use the adv. DPF, 13 out of 26 IMBH binaries in the galactic 
globular clusters lie withing the observable reach of the detector. This makes the 
event rate larger than Eq. (9) by a factor 13/2 — 6.5. 

For the latter case, currently, more than 10 galactic massive young clusters 
(GMYCs) have been discovered^^^. Giirkan et al}'^'^ performed numerical simula- 
tions and found that two IMBHs may form in GMYCs if the initial binary fraction 
is relatively large. After IMBH binary formation, it shrinks due to the dynamical 
friction with the cluster stars. The timescale of the dynamical encounters would 
be < IGyr For simplicity, we assume that IMBH binaries are all located at 
the galactic center. From Fig. 5, one sees that DPF is not sensitive enough to de- 
tect GW signals from IMBH binaries in GMYCs. On the other hand, when we use 
adv. DPF, it has ability to detect all of the IMBH binaries in GMYCs. Following 
Fregeau et al.^^^, we assume that the number of star clusters massive enough to 
form IMBH binaries is the same as the one of globular clusters, and 10% of them 
actually produce IMBH binaries. In our galaxy, about 150 globular clusters^^""^ have 
been observed. This means that at least 15 IMBH binaries are expected to be within 
the reach of adv. DPF. The detection rate of the IMBH binaries in GMYCs with 
adv. DPF can be roughly estimated as iV^MYC ~ 15/(13. 8Gyr) = 1.1 x lO-^yr"!. 



2.4.2. IMRIs 

After supernova explosions, many stellar-mass BHs of ^ IOMq would form. These 
would sink to the cluster cores due to the mass segregation, and their number 
density is expected to be comparable to the one of the main sequence stars, typically 
n = lO^pc"'^. In this subsection, we consider IMBH binaries with intermediate-mass 
ratio inspirals (IMRIs). In Fig. 8, we show observable ranges of DPF (red thick 
solid) and adv. DPF (green thin solid) with 7712 = IOMq. One sees that DPF is not 
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Fig. 8. Observable ranges against total mass M with m2 = IOMq . The meaning of each 
curve and plot is same as Fig. 5. 



sensitive enough to detect any GWs from IMRIs with TO2 = IOM0, but adv. DPF 
has abihty to detect these signals from some of GCs and the galactic center. 

For GMYCs, from Fig. 8, one sees that even adv. DPF is not sensitive, on aver- 
age, to an IMRI signal of (10 + 1O^)M0 at the galactic center. We found that 1.7% 
of these binaries have SNRs greater than 5 due to the optimal orientations^^^. On 
the other hand, the merger rate due to three-body interaction per cluster becomes 
^GMYC ^ 2 X 10"*yr-^ 213 jj^g^j^g conservative assumption that 50 GMYCs 

may contain IMBHs (see e.g. Refs.^^^'^"'^^). In this case, the total detection rate 
becomes 

^GMYCtot ^ Q Q^j X 50 X 2 X IQ-V"^ = 1-7 X 10~V~^- (10) 

For GCs, we found that the merger rate is limited by the requirement that the 
merger rate in a single cluster cannot exceed Nmass = ('^i/'^2)T'ago ^^^'2'^®, where 
Tage is the current age of the cluster. Otherwise, the mass of the larger BH has been 
increased considerably. From Fig. 8, one sees that IMRI signals from about 3 out 
of 26 globular clusters might be detected by adv. DPF. This leads to the estimate 
that it can detect IMRI signals from about 3 x 150/26 = 17.3 globular clusters in 
total, and hence the detection rate of globular clusters in total becomes 

^GC.tot < ^ 2 X 10"V~^- (11) 



■^For relatively large IMBH masses, the two-body interaction may give comparable contribution 
compared to the three-body interactions^*^®' . 
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3. DECIGO/BBO 
3.1. Basic Designs 

DECIGO is first proposed by Seto, Kawamura and Nakamura^^. It consists of four 
triangular sets of detectors whose configuration is shown in Fig. 1. Tliis effectively 
corresponds to eight individual interferometers. Each triangular detector has an 
arm-length of lO'^km. Its primary goal is to detect PGWB with the GW energy 
density of ficw = 10~^^. Since the WD/WD confusion noise will have a cutoff 
frequency at around 0.2 Hz ''^ . DECIGO has an advantage on detecting this source 
over LISA (see Refs.®^'^^'^^'* for the discussions of the NS/NS confusion noise). 
Therefore two of the four triangular detectors are located on the same site forming 
a star of David so that correlation analysis^-'^^^ can be performed to detect PGWB. 
The rest of the two detectors are placed far apart to increase the angular resolutions 
of the source locations. 

BBC has almost the same constellation as DECIGO. The main difference is that 
while DECIGO is a Fabry-Perot type interferometer, BBO is a transponder-type 
interferometer. BBO has arm-lengths of 5 x lO^'km. 



3.2. Noise Spectrum 

The noise spectrum of BBO is given as follows. The non-sky-averaged instrumental 
noise spectral density for BBO is obtained from Ref.^* as®^ 



sr\f) 



1-49 



1.8 X 10""'-' ( —\ + 2.9 X 10""^^ + 9.2 x 10"^^ ( — ^ 



-4 



Hz-1 



iHzy viHzy 

(12) 

It has 20/3 times better sensitivity than the one for the sky- averaged sensitivity 
Apart from instrumental noise, there are astrophysical foreground confusion noises. 
These confusion noise spectral densities and the energy densities of gravitational 
waves are related as^'' 

sr' = §^r'PcnGw (13) 

Here, pc = is the critical energy density of the Universe and 

_ 1 dpGw ,x 

"GW^— -7j — T (14) 

is the energy density of GWs per log frequency normalized by pc. The energy den- 
sity of GWs that originate from extra-galactic WD binaries has been estimated as 
17gw ~ 3.6 X 10~^^(//10^''Hz)^/'^ which leads to the noise spectral density of^ 

S!^^-s^\f) = 4.2 X 10-^' ' Hz-i. (15) 
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On the other hand, the one from galactic WD binaries has been calculated as^^ 

-7/3 



ST\f) = 2.1 X lO-'i^ (^jL.^ Hz-i. 



(16) 



We muhiply 5,^ ' and S'^^^s'^i by a factor T = exp{-2 (//0.05Hz)^}, which corre- 
sponds to the high frequency cutoff for the white dwarf confusion noises. We also 
have to take into account the confusion noise from NS binaries. Its noise spectral 
density S^^ is estimated as^^'^-^^ 



5r(/) = 1-3 X 10- 



1 Hz 



-7/3 



no 



10^7 Mpc^Vr^^ 



Hz" 



(17) 



where no denotes current merger rate density of NS/NS binaries. Putting altogether, 
the total noise spectral density for BBO becomes 



Snif) = 



Sir' if) 



cM-^T-^ldN/dfY 

+ 5r"8^H/)-^(/)+-^clcan5r(/) 



if) + ST\f)T{f) 



(18) 



Here dN/df is the number density of galactic white dwarf binaries per unit fre- 
quency, which is given by^^'^ 



dN 



= 2 X lO^^Hz"^ 



/ 



1 Hz 



-11/3 



(19) 



K ~ 4.5 is the average number of frequency bins that are lost when each galactic 
binary is fitted out. The factor J-c\can in front of S^^{f) represents our assumption 
of the fraction of GWs that cannot be removed after foreground subtraction. In 
this review, we assume that NS/NS foregrounds can be subtracted down to the 
level below the instrumental sensitivity. The lower and higher frequency ends of the 
BBO sensitivity band are set as /low = 10~'^Hz and /high = lOOHz, respectively. 
The noise spectrum of BBO is shown as a red thick solid curve in Fig. 2. 

DECIGO has been proposed with 3-4 times less sensitive spectrum than BBO. 
Its instrumental noise spectrum is given by 



iJ DECIGO (/) = 5.3 X 10 



-48 



(l+2/') + 



2.3 X 10^^ 2.6 X 10" 



Hz- 



(20) 



where y = ,f / fp with fp = 7.36Hz. However, this is not the fixed design sensitivity 
and there is a project going on to improve the sensitivity to the same level as BBO. 
The frequency range of DECIGO/BBO is given as 



(/low I /high) 



(10"^Hz,10^Hz) 



(21) 
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3.3. Event Rate 

The promising sources for DECIGO/BBO are inspirals of NS binaries. In this sub- 
section, fohowing Cutler and Harms^'', we derive the detection rate of these sources. 
Here, we adopt the model where NS/NS binaries only exist at redshift below z = 5. 
(This is inferred from observations. See the discussion below.) Since DECIGO/BBO 
has enough sensitivity to detect the farthest NS/NS binaries considered here, the 
merger rate corresponds to the detection rate. The NS/NS merger rate can be writ- 
ten as 

poo J 

iVNSNS= / ^TT[aori{z)fn{z)-^dz, (22) 



where 



, , /"^ dz' dr 1 

Jq H[z) dz (l + z)H(z) 

Here, H{z) represents the Hubble parameter at redshift z. For A-CDM cosmology, 
it is given by 



H{z) = Ho^{l + z)^n^ + flA. (24) 

n{z) is the merger rate of NS/NS binaries at redshift z and it can be re-expressed 
in terms of the current merger rate no and the redshift dependence R{z) as n{z) = 
ho X R{z). For R{z), we adopt the following piecewise linear fit ^'^'^is based on 
observations discussed in Rcf.^^^; 

(l + 2z (z < 1) 
R{z) = i |(5-z) (1 < z < 5) (25) 
lo (z>5). 

Then, we finally obtain the detection rate as 



VVnsns = 10' ( "„3 ) yr^i . (26) 



,10-'^Mpc~Vi-\ 

The detection rate for BH/NS binaries are expected to be roughly 1/10 of NS/NS 
detection rate. 

We can also estimate the NS/NS foreground S^^{f) as follows. First, we apply 
the convenient formula given by Phinney^^'^ as 

9 H^,^ f J, ^^1 + Z)4/3H(Z) 

. 3.74 X 10-^.,3 (My- n)-^ r,,-^) . (27) 



By using Eqs. (13) and (27), we obtain 
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Fig. 9. A circular binary with masses mi and m2 . 



4. Gravitational Waves from Compact Binaries 

In this section, we focus on the GWs from inspiral binaries composed of compact 
stars sucli as BHs and NSs within the context of GR. In Sec. 4.1, we first derive the 
GWs from a circular binary at the leading order. Then, In Sec. 4.2, we extend the 
result to higher order terms. Finally, in Sec. 4.3, we consider binaries with small 
eccentricity. 



4.1. Quadrupole Gravitational Radiation from Binaries with 
Circular Orbit 

In this section, we derive the leading order gravitational waves from compact bi- 
naries, namely the quadrupole radiation. The binary circularizes quickly during its 
orbital decay due to the gravitational radiation reaction, and hence we restrict our 
attention to the binary with circular orbit in this section. Also, we assume that the 
velocity of each body is sufficiently small compared to the speed of light. 

Let us assume that we have two compact objects with masses mi and m2 in the 
x-y plane with orbital separation a at a distance from the detector. We define 
the total and the reduced mass as Mt = mi + m2 and ^ = mim2/Mt, respectively. 
The binary components obey the Kepler's Law at the leading order with orbital 
angular velocity defined as 




As shown in Fig. 9, we denote the positions of these objects in the center of mass 
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frame as (xi,yi) and (x2,2/2), which can be expressed as 

(a;i,yi) = (^acos(r!i+|) ,^asin(m+|)^ ' ^^^^ 

By taking the time derivative of the quadrupole moment Mij — X)a=i 2 '"^aXai^aj 
twice, we obtain 

h+ = A+cos2Vtt, (32) 
hy, ^ Ay, sin2r2t . (33) 

Here, we have defined 

A^^'j^iL-n). (35) 

where i = £ is the unit orbital angular momentum vector and n is the unit vector 
pointing towards the direction of GW propagation. From this result, we see that 
the monochromatic GW frequency can be written as / = ri/vr at the leading order. 

Next, we take the effect of radiation reaction into account and estimate the GW 
phase 4){t) = / 2nfdt First, the total energy of the binary system is given as 

Next, by using the quadrupole formula 

^ = ^(M.,M.-i(M..)^), (37) 

we have the radiated energy due to GWs as 

dEG^ 32 iP-Ml 



dt 5 a5 ■ 

By using the balancing equation 

dE dEoYf 
'dt ~ df ' 
the evolution of the binary separation can be expressed as 



(38) 



(39) 



a dE 64 riM? , , 



which can be solved to yield 



a(t)^[—i,Ml(t^-t)\ , (41) 
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where to is the coalescence time. This can be turned into the frequency evolution 
as 

I _^ _ 3 M/^^ a 
TT 2 TT 

= ^^8/3^5/3^11/3^ (42) 

5 

where A4 = Alf^^n^^^ is the chirp mass. This equation can be solved as 

/ 5 \ 3/8 I I 
■^^\256) nM^/^ {to - ' ^^^^ 
By integrating this by once, we get the phase as 

q^it) = J 27Tfdt = -2 Qx-i(to -t)y + 00, (44) 

with 00 representing the coalescence phase. 

Next, we derive the waveform in the Fourier domain h{f) defined as 

/oo 
e^'''f%{t)dt. (45) 
-co 

Since we have assumed that the velocity of each binary component is much smaller 
than the speed of light, the amplitude A{t) and the phase satisfy the following 
conditions: d\nA/dt <C d(j)/dt and d'^<j)/dt^ <C {dip/dtY. Under this situation, we 
can apply the stationary phase approximation, where the Fourier component of a 
function B{t) = A{t) cos (p{t) becomes^^ 

B{f ) « Uit) (i) cxpj(27r/t - 0(/) - n/A). (46) 



2 \dt_ 

Here, t is the time that satisfies d(f>{t)/dt = 27r/ and 0(/) — 0[i(/)]. From Eqs. (43) 
and (44), we have 

t{f)=to-^M{STiMfr^'\ (47) 
0(/) = 00 - 2(8^X/)-5/3 . (48) 

Therefore, we obtain 

K{f) = (^i±i^ j e^*+(/), (49) 

/ix(/) = ^r'/'(i-n)e**x(^), (50) 
where the amplitude and phases are given as 

IT 1 

*+(/) = 2^/to-0o- J + *opn(/), (52) 
vI/x(/) = *+(/) + |, (53) 
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respectively. Here, ^opn(/) is defined as 



and Dl is the luminosity distance given as 



Dl^{1 + z) I 
Jo 



dz' 



(54) 
(55) 



4.2. Higher Post-Newtonian Corrections 

In the previous section, we derived the waveform at the leading order, where wc used 
the quadrupole formula for the radiated energy flux and assumed that the binary 
orbit is of Newtonian. In this section, we extend the previous analysis to higher 
orders. Since we go beyond Newtonian, this expansion is called the post-Newtonian 
(PN). Our expansion parameter is (tj^/c^) where v is the typical velocity of the 
binary component, and the terms that are proportional to (u^/c^)" are of nPN 
orders. In this section, we consider the extension up to 2PN order^^^. 

Furthermore, we neglect the PN corrections to the amplitude. This is because 
when we take the correlation between two different waveforms, it is more sensitive 
to the deviation in the phase than the one in the amplitude. This type of waveform 
is called the restricted PN waveform. 

In Rcf.^^^, expressions for E, E and 51 are shown to 2PN order. By using them, 
up to this PN order, / becomes 



/ ^ ^7r8/37W5/3yii/3 
5 

/ 34103 13661 



1 - 



V 18144 



743 
336 



59 2 

77 H 77 

2016 ' 18 ' 



^77 ) x + (4^ - /3)xi/2 



(56) 



where 77 s fi/Mt is the symmetric mass ratio with PN expansion parameter x e 
(TrMt/)^/^ — Alt /a — ■ The spin-orbit and the spin-spin couplings are given as 

^2 



12 ^ 

1=1,2 



113- 



7bri)L- Xi, 



^.{-247x. 



X2 + 721 [L-xi] (L-X2 



0} 



(57) 



(58) 



where the dimensionless spin parameter is defined as Xi = Si/mf with Si denoting 
the spin angular momentum of the i-th compact object. The gravitational wave 
phase in the Fourier domain can be calculated as 

/3715 



*+(/) = 27r/to - 00 - ^ + *opn(/) 
-4(47r - /3)a;^/2 ■ 



1 



15293365 



508032 



V 756 
27145 
504 



55 , 



-11 



3085 
^72~ 



7^2 - lOcr ] 



*x(/) = *+(/) + 2 ■ 



(59) 
(60) 
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4.3. Gravitational Waves from Binaries with Slightly Eccentric 
Orbit 

In this section, we consider gravitational waves from a binary that have a small 
eccentricity e. We can express the eccentric orbit as 

1 + ecosi2t 

where the orbital angular velocity can be written as 



n ^ VWEZl. (62) 



By using the quadrupole formula, the radiated energy flux can be estimated as 

.2 Ji,r5 



dt 5 a5(l-e2)7/2 1^ +24'' 96"" 
and similarly, the radiated angular momentum flux is calculated as 



(63) 



— ^^ijk \Q xyiQyy Qxx) ^QxyiQxx Qyy) 

Next, we derive the evolution of the semi-major axis a and the eccentricity e. 
Since the orbital energy of the binary is given in Eq. (36) and the angular momentum 
is given as 



2 



J, = ^ri^Q = ^ Va(l-e2) (55) 

at the Newtonian order, by using Eqs. (63) and (64), the evolutions of a and e can 
be calculated as 

da da dE 
'dt ~ 'dE~di 

64 r^Mf / 73 2 37 4 



5 a3(l-e2)V2 y 24 96 

de 1 — f 1 da 1 dJ^ 
- 2- 



2e \a rfi 

304„ .M.^ ^..121^2). (6,) 



15 a4(l-e2)5/2 y 304 
From these equations, we obtain 



dn 19 1 + — + ^e* 

da _ 19 1+24^ +96^ (gg) 



« 12e(l-e2)(l + i§e2) 



Therefore, by integrating once, we obtain the relation between a and e as 

(69) 



a _ 1 - ef 


/ \ 12/19 
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Oi 1 — 




.u 
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304 i J 
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where and denote the initial semi-major axis and eccentricity, respectively. 

Finally, we derive the gravitational waveform (in the Fourier domain) from the 
compact object with the small eccentric orbit^^'^^^. By substituting Eq. (69) into 
the GW frequency / = JI/tt, we obtain 

/(e) = /.^, (70) 
C(e) 



ifip^A 22^ 

(1 - e2)3/2 

where fi is the frequency at e = e^. Especially, when e ^ 1, we have 

/ .X -19/18 

so that up to 0{e), the asymptotic eccentricity invariant, 

le = e^r"" (73) 

is conserved. 

The evolution the frequency / is given by 

/ = ^TT^/^A^^/^"/^ (l + ^e^) + 0{e^) , (74) 

for e <C 1. By substituting Eq. (72) into the above equation, we have 
96 



^ 5 



,^8/3_^5/3^^11/3 ^ ^J^J-19/9^ _ (75) 

From this equation, we can calculate t{f) and </>(/), leading to the GW phases in 
the Fourier domain as 



*+(/) = 2^/t(/)-0(/) ^ 



2^ ft. - 00 - 5 + *OPn(/) (1 - ll^e/-^^/^) 



(76) 



vl/x(/) = «-+(/) + !. (77) 



5. Observed PN Inspiral Waveforms 

In this section, we explain the response of two-armed interferometers to GWs from 
compact binaries within the context of GR. Expressions below can be applied to 
LISA and DECIGO/BBO. (See Sec. 2.3 for DPF.) As explained in the previous 
section, we use the so-called restricted PN waveforms where we keep the higher 
PN terms in the phase but only take the leading terms (leading PN and dominant 
harmonic contributions) for the amplitude^^. In Sec. 5.1, we consider the case where 
the orbital angular momentum L and the spin vectors Si {i ~ 1,2) are all aligned 
(or anti-aligncd). Then, in Sec. 5.2, we include the effect of precession and see how 
this modifies the waveforms. Finally in Sec. 5.3, we explain the fitted waveforms 
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including all of inspiral, merger and ringdown phases for the spin-aligned binaries 
with circular orbits. 



5.1. Spin- Aligned Binaries 




Fig. 10. We use two types of coordinates: (i)a barred barycentric frame (x, y, z) tied to the ecliptic 
and centered in the solar system barycenter, and (ii) an unbarred detector frame (x, y, z), centered 
in the barycenter of the triangle and attached to the detector. (This figure is taken from Ref.-"^^".) 



Following Ref.^^^, we introduce two Cartesian reference frames: (i) a barred 
barycentric frame (i, y, z) tied to the ecliptic and centered in the solar system 
barycenter, with z (unit vector in z direction) normal to the ecliptic and xy-plane 
aligned with the ecliptic and (ii) an unbarred detector frame {x,y,z), centered in 
the barycenter of the triangle and attached to the detector, with z normal to the 
detector plane (see Fig. 10). The orbit of the detector barycenter in the barred 
frame can be written as 

e{t) = 7r/2, 4){t) ^ 27rt/r, (78) 

where T = 1 yr and we have assumed (j)(t = 0) ~ 0. 

A triangular detector has effectively two individual interferometers. Therefore, 
it is possible to measure both polarizations with a single detector. We first focus on 
the interferometer I which consists of the arms 1 and 2. This detector measures 
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where SLi{t) and 5L2{t) are the differences in lengths in arms 1 and 2 due to GWs. 
L is the arnilength where GWs are not present. The factor -\/3/2 conies from the 
fact that the opening angle of adjacent arms is 60°. Next, we introduce two principal 
axes for the wave; p = N x L/\N x L\ and q = p x TV, where L is the unit vector 
parallel to the orbital angular momentum and N is the unit vector pointing towards 
the center of mass of the binary. The two polarizations are exactly 7r/2 out of phase 
and the waveform becomes 

habit) = A+H+coa^{t) + A^H^^sin^it), (80) 

where and H^f^ are the polarization basis tensors, defined as 

^ab = PaPb - qaQb, = PaQb + QaPb- (81) 

From Eqs. (79)-(81), the detector output hi{t) becomes 

fn AT 

hlit) = ^A+F+iOs, 0s, V-s) cos0(i) + ^A^F,^ {0s,^s,i^s) sin0(i). (82) 

Here, Fj^(0s7 ^s, "^s) and F^ {6s,4's,4's) are the detector beam-pattern coefficients 
for the plus and cross polarization modes, respectively, and when the detector is an 
interferometer, they are given by 



F+{9s,(f>s,i^s) = + cos^ 0s) cos(20s) cos(2V's) - 008(6*3) sin(2(/.s) sin(2V's), 



2 



(83) 



i^i^ (6(3, 0s, V's) = J(l + cos^ ^s) cos(20s) sin(2V^3) + cos(6ls) sin(20s) cos(2V;s)- 



2 



(84) 



{Os-,4>s) represents the direction of the source in the detector frame and '0s is the 
polarization angle defined as 

tan^s = ^^-^-^-^^:^^^^-^^ (85) 
p-z N-{Lxz) 

Also, another interferometer can be constructed from arms 2 and 3. We call 
this interferometer IF and its signal can be written as hw = {SL2{t) — 6L3{t))/L. 
However, since hi and hn' have some correlations, they are not independent in- 
terferometers. We combine interferometers I and H' to construct interferometer H 
which is uncorrelated with detector L The signal of interferometer H is 



huit)^^Mt)+2hn'it)]^^ 



„ {hxy + hyx) 



(86) 



This interferometer H corresponds to the one that is rotated by 45° with respect to 
interferometer L Thus the beam-pattern coefficients for the interferometer II are 

^il (^s , 0s , ^s ) = (^?s , 0s - 7r/4, 0-s) , (87) 
^i^ (^s, 0s, V-s) - F,"" {0s,(f>s - ^/4, ^s). 
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When wc perform parameter estimation, we include the direction (0s, 0s) and 
the orientation ((?l7'/'l) of the source measured in the solar barycentric frame in 
binary parameters. Therefore we need to express the angles Os{t)^ 4's{t) a-nd ips{t) 
in terms of 0s,(j)s^ Oh and 0l- Explicit formulas are given in Appendix A.l. 

Re-expressing the waveform measured by each interferometer in terms of an 
amplitude and phase, it becomes 

ha (t) = ^ Apol.a (<) COS + V3pol,a (t) + fD (0] , (89) 

where a = I, II is the label to specify the interferometer. The polarization amplitude 
^poi,c<(i), the polarization phases (p-po\,a{t) and the Doppler phase ifioit) are defined 
as 



Sin((ppol,a(0) 



{l + {L-Nf)F+{t) ^ 

^pol.Q it) 

2{L-N)F-{t) ^ 



ifoit) - 27r/(t)i?sin0scos[<^(i) - 0s 



(90) 
(91) 

(92) 
(93) 



where R = lAU. The Doppler phase denotes the difference between the phase of 
the wavefront at the detector and the phase of the wavefront at the solar system 
barycenter. It arises from the fact that the detector orbits the Sun. 

Later, we estimate the measurement accuracies of the binary parameters using 
the matched filtering analysis, where we work in the Fourier domain. Therefore we 
calculate the Fourier transform of the signal using the stationary phase approxima- 
tion to yield^^^ 



^GR 



(/) 



V3 



Af- 



-7/6gi*(/) 



r^pol,a(i(/)) 



-«(Vpo1,o+'(=d) 



where the amplitude A and the phase 5'(/) are given by 

1 M^/^ 

A = 



30^2/3 Dl 



(94) 



(95) 



*(/)^2./to-0o-5 + 4(.A^/)-/^[l-^ 

/3715 55 \ ^.^ ^, 
+ ( ^ + —vj X ~ 4(4^ - /3)x3/2 

/ 15293365 27145 3085 , „ \ 2 



V 508032 



504 



72 



(96) 
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Also, when wc take the average of the waveform (Eq. (94)) over the angles, the 
sky-averaged (or pattern-averaged) waveform becomes 

^GR(/)-^^r'/^e'*(^^. (97) 



2 



5.2. Precessing Binaries 



In this section, we introduce an additional effect, the precession. The spin-orbit 
interaction and the spin-spin interaction change the orientations of the orbital an- 
gular momentum vector L and the spin vectors Si. These vectors precess over a 
time scale longer than the orbital period but shorter than the observation period. 
This effect drastically changes the detected waveforms. 

5.2.1. Simple Precession 

The precession equations are given in Appendix A. 2. In this review, we assume that 
one of the spins of the binary constituents is negligible (i.e. S'l ~ 0). Then, there 
do not exist spin-spin interactions. We also assume that the orbital angular mo- 
mentum L is neither parallel nor anti-parallel to the total spin angular momentum 
S {= Si + 82)- Then, the precession equations are simplified and L is obtained 
analytically up to some approximate orders. This is the so-called simple precession 
approximation^'^'^ . (This also holds when the masses of the binary constituents are 
equal (mi ~ and spin-spin interactions are negligible, instead of Si ^ 0.) 
Under this approximation, the precession equations simplify to 



L=[2 + :,^]-,xL, (100) 



i = -^- — , (98) 
b a \ a J 

5 = 0, (99) 
3 m2 \ J 
2 mi J a'^ 

^-(2 + |!!^)4x5, (101) 

where J is the total angular momentum J = L + S. Under this approximation, the 
following quantities become constant during the inspiral. Si ■ S2, n = L ■ S, and 
the magnitude of the total spin angular momentum 5* = jS'i -I- S'2|. 
Next, we define a quantity "f{t) as 

7(0 - (102) 

Then, the magnitude and the direction of the total angular momentum J can be 
expressed in terms of k, 7(t), L{t), L and S as 



J = LVT+2^^7 + 7^, (103) 
V 1 + 2k7 + 72 
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From Eqs. (100) and (101), it can be seen that both L and S precess around J 
with the precession angular velocity 

a,^(2 + l^-^)±,. (105) 
\ 2 mi / a"^ 

In general, the processing time scale ilp^ is shorter than the radiation reaction time 
scale L/\L\. Therefore from J = LL, J changes in magnitude but the direction is 
almost constant. Actually, if J is much smaller than L (as this can happen when L 
and S are anti-aligned with almost the same magnitudes), J can change significantly 
in one precession period. Therefore, we introduce the following small parameter. 

Then, J processes around the fixed direction Jq with 

Jo = J -eJ X L + 0{e^). (107) 
To the same order, the precession equation for L becomes 

L = QpJo X L + eflpiJo x L) x L + 0{e^). (108) 

The solution of this equation can be obtained geometrically^^''. Let us assume that 
Jo points in the (^j , ) direction. We denote Al as the opening angle of the cone 
on which L processes (i.e. the angle between L and Jq; sec Fig. 11). This can be 
regarded as the angle between L and J apart from the errors of order and is 
given by 

cosAL=i-J= , ^^'V - , (109) 

VI + 2k7 + 72 



\L\ jVT^ 



2 



sin Al = ^ = , ' \ (110) 
^1 + 2^7 + 72 

Then, L can be expressed as 

T T \ , Z-Jocos9] -J ^ sinALsina 

L = Jq cos Al H : — ^ sni Al cos a + (Jo x z) — = , (HI) 

sin 9 J sin Wj 

where a is the precession angle defined as the solution of 

We assume that a = is realized when i • 2; is maximum (see Fig. 11). By solving 
the above equation, we obtain 

- 5 1 I ^"^1 

2(gL)^ - 3kS(L + KS)gL - 3kS^(1 - K2)arcsinh ( — -— - 

\SVl - 




Fig. 11. The unit orbital angular momentum vector L precesses around a fixed vector Ja ^"^^ . 
The opening angle of the precession cone is given by Al and the precession angle is denoted as a. 

where Uc is a quantity which characterizes a at t — tc and Q is defined as 

g = v/1 + 2k7 + 7^- (114) 

5.2.2. Detector Response 

When precession is taken into account, the principal axis p varies with time so that 
the GW phase (j>{t) no longer equals to twice the orbital phase 4>oTh{t) = / fl{t)dt. 
We define this diff'erence 2S(j){t) by 

^{t) = 2(j>o,b{t) + 26(j){t). (115) 

This S(l){t) is the so-called Thomas precession phase. We absorb the constant of 
integration in the definition of (l)oib{t) so that 4>oih{tc) = (f'c- In general, S(j){tc) ^ 
so that (f)c docs not equal to 0(ic) anymore in this case. From Eqs. (82) and (115), 
the detector output ha{t) becomes 

ho.it) = ^A+F+{es,(bs,^s)cos2{<j)orhit) + 5(b{t)) 

+^A^F^ {0s,4's,i^s) sin2(</)o,b(i) + S^it)) 

= ^^4?(^"'co«20orb + Fr sin2</)o,b), (116) 
2 aD 
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where F^°^ and are defined as 

F^°' = {1 + {L- Nf)F+ cos25(j) - 2{L ■ N)F^ sin 2(5^, (117) 

F^'" = -{1 + {L- Nf)F+sm264>-2{L- N)F^ cos2S(t). (118) 

Following Eq. (89), we express this output (116) in terms of an amplitude and phase 
form, and we also take the motions of the deteetors into account to obtain 



K{t) 



2 aD 

where A^^'^^''^^{t) and ^I'^ai^) gi^^n by 



cos 2^,,^{t) + ^PXa(i) + VDit) 



Frit) 



Frit) 



The Fourier transform of this waveform is given by 



jK':lw)) 



-''■{Vl'al,^+VD) 



(119) 

(120) 
(121) 

(122) 
(123) 



which looks similar to Eq. (94), except Apo\,a{t{f)) and ippo\,a{t{f)) are replaced 
with K'ohaW)) and V^l^oLW)), respectively. 

Finally, we need to calculate the Thomas precession phase S(j){t). Apostolatos et 
al?'^'^ derived an explicit form as 



m) 



dt 



L-N 



1 



{L xN)-L. 



(124) 



(i- Af)2^ 

This expression includes an integration which makes the computational time very 
long. Vecchio^^^ estimated the binary parameter accuracies by choosing a few ran- 
dom sources with and without including (50(i) and found that this term did not 
affect the results, and concluded that it could be neglected. This is true for the 
binaries for which L ■ N never becomes close to ±1. However, when i • TV « ±1, 
the direction of the principal axis p changes rapidly with time, and hence the po- 
larization angle V's {t) also changes rapidly. Since it is the Thomas precession phase 
5(j}{t) that cancels this rapid change, 54'{t) cannot be neglected in this case. 

The direction which 5(j) is measured from is arbitrary. It can be seen from 
Eq. (124) that Apostolatos et al.^^'^ defined Sep as the angle measured from the 
principal axis at the time of coalescence. In general, since p\t=to does not lie in 
the orbital plane at a time t, we need to follow the evolution of the principal axis 
to calculate the Thomas precession 5<f>{t). This is the reason why the integration 
appears in Eq. (124). 

Here, we try to derive an approximate expression of 54>{t) that is not in the 
integral form. To do so, we use a specific vector L x ii that always lies in the orbital 
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Fig. 12. The Thomas precession phase &(t>{t) is the angle from the vector ii to the principal axis 
p = N X L ^''^'^ . Wo define 5(f)' as the angle from the vector i X li to the one p. Notice that these 
two vectors always lie in the orbital plane. S(j>' equals to Sep + 7r/2 up to 0(\SL\/L). 



plane, where tit is a constant unit vector. In the Hmit of a large separation (i.e. 

ooorL— >cxj),i approaches to Jq. Using this Jq, we take u = N x Jq 
(see Fig. 12). We define an approximate Thomas precession phase S(j){t) by the 
angle from this vector ii to the principal axis p. (We may take the vector u to be 
the principal axis at the time of coalescence p\t=to- Then, (50(io) becomes and 
4>{to) — (pQ. However, the computation is easier for choosing it to be iV x Jq.) We 
denote the angle from i x u. to p as Scfi' , and the difference between Jq and L as SL; 
L = Jo + 5L. Then, Scj)' can be related to 64> as 64>' = 64> + tt/2 up to 0{\6L\/L). As 
64>' always lie in the orbital plane, by using this angle we can express S(t>(t) without 
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integration, sin 50 and cos (50 can be written as 
sin (50 « — cos (50' 



(125) 



{Lx N)-{Lx u) 
\L X N\\L X u\ 
{L-N){L-u) 

(l-iL-Nmi-iL-ur) 

^5(1)' 

(126) 

'(i-(i,-iv)2)(i-(i,.{t)2) 

Therefore the Thomas precession phase d(j){t) can be expressed as 

50^arctanf-(-^:^j^-^-")V (127) 
Here, the exphcit forms of L ■ N, L ■ ii and N ■ {L x ii) are given in Appendix A. 2. 
5.3. Inspiral- Merger- Ring down Hybrid Waveform 

In the previous 2 sections, we explained how inspiral waveforms can be modeled 
using PN theory. One possible extension would be to include merger and ringdown. 
The phenomenological inspiral-merger-ringdown hybrid waveforms have developed 
by Ajith et aZ.^^^'^os. They first performed numerical simulations to generate merg- 
ing binary BH waveforms from 8 cycles before merger including non-precessing BH 
spins. Then they matched these waveforms with the Taylor Tl post-Newtonian 
(PN) inspiral waveforms^^^ at 3.5PN phase accuracy^^''. They included the 3PN 
amplitude corrections^^^^^ao g^^^^j -(.j^g 2.5PN corrections coming from the spins^'^^ to 
the dominant quadrupole mode (see Ref.^^^ for detailed procedure of the matching). 
The Fourier component of the (observed) matched waveform 

h{f) ^ yl(/)e'{*(^)-'^p°'(^)-'^°(-^)> (128) 

is fitted by using phenomenological parameters as 

[ /'-^/^(i + El2«.^'/') (/</i) 

A{f) ^ C{t{f))f;"' + j^l^ (/^ </</,), (129) 

[wrCif,h,a) (/2</</3) 



and 



«-(/) = 27r/to - 00 - ^ + *0PN ( 1 + II ^''^"'/'fe ) (130) 



7 



k=2 



"^The polarization phase V'poiC/) and the Doppler phase foif) are given in Eqs. (91)-(93) without 
the index a. 
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with /' = ///i- /i and /2 are the transition frequencies from inspiral to merger 
and merger to ringdown phases, respectively, with Wm and Wr corresponding to the 
normaUzation factors to make A{f) continuous at these frequencies, corresponds 
to the cutoff frequency. C is a numerical factor discussed in the next paragraph, 
are given as 

ei = 1.4547X- 1.8897, (131) 
£2 = -1.8153X+ 1.6557. (132) 

Here, the spin parameter x is defined as 

X = Xi + — 2— X2, (133) 

where 5 = {mi — 7712) /A/j and Xi ^ Si/ra^ with and Si being the mass and 
the spin angular momentum of i-ih BH, respectively. For the equal spin case, x 
becomes X = Xi = X2- On the other hand, are obtained from (Z,m) = (2, ±2) 
FN waveform^^^ as 

323 451 , , 

"^ = -224 + 168^' (1^4) 

/27 11 \ , , 

a3^(^-^.jx. (135) 

£(/, /2,o') is the Lorentzian function centered around the frequency /2 with width 
a. The phenomenological parameters -0^ and = (/ii /2j cr, /a) are fitted against 
physical parameters {Mt,r],x) as 

3 N 

^'^-EE4"''^^v+^.°, (136) 

i=l j=0 

t'^'l^l^ ' (137) 

i=l j=0 

where N = min(3 — i, 2) while Vk ''^ ^ V^fe a-nd /i° are listed in Table I of Ref.^*^^. 

The numerical factor C{t{f)) is given by^^^ 



where Apo\{t{f)) is given in Eq. (90) without the index a. For the time t(f), we 
take up to 3.5PN order. The expression for df / dt up to this order is given in Arun 
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et al.' 



Wc integrate this to yield 
5 



tif) = to 



256 



4 /743 
3 V336 



11 , 

+ —7] X 



5 
Stt 

Y 



(47r - I3)x^l'^ + 2 



/ 3058673 5429 617 



V 1016064 



77 — a \ X 



7729 



672 

/ 15335597827 



^r^]x'/' + l5[- 



616A 45l7r2 



1008 ' 144 
10817850546611 



V 60963840 



25565 



-V 



17127 327r2 



5184 ' 105 3 
154193357r 757037r 



48 
3424 
^ 1575 
148097r 



93884313600 

15211 



1016064 



6048 



-77 



3024 



6912 
log (32768a;)) 
7/2 



\x- 



(139) 



Here, 7 = 0.577 is the Euler's constant, A = and 9 = We take the 

contribution of spins into account up to 2PN order with the spin-orbit coupling /3 
and the spin-spin coupling a. The sky-averaged values of the beam-pattern functions 
are (F^^) = (F^^) = 4/15 which yield the sky-averaged value of C as (C) = 

We use this hybrid waveform when we perform binary parameter 



estimation with DPF in Sec. 8. 



6. Brans-Dicke Theory 

Currently, solar system experiments put stringent constraints on gravitational theo- 
ries in the weak field regime^^^, but general relativity (GR) is not the only one that 
passes this test. Among the modified gravity theories that pass these tests, in this 
section, we consider one of the simplest extension of GR called the scalar-tensor thc- 
Qj-yi59^ As its name tells, in this theory, we add scalar degrees of freedom to gravity. 
Among various types of scalar-tensor theories, we here focus on Brans-Dicke the- 
ory^^*. Wc note here that the scalar-tensor theory had already been proposed 20 
years before the paper by Brans and Dicke^'^'^'^'^*. 

The reason we consider scalar-tensor theory is that it may give clues to solving 
unknown problems in GR and cosmology: (I) Current observations show the acceler- 
ation expansion of the Univcrse^^^. However, the cause of this cosmic accelerating is 
unknown. If we are to explain this phenomenon with a matter field, we need to intro- 
duce the dark energy which has a negative pressure^^°'^^^. Another approach is to 
modify gravitational theory from GR. In scalar-tensor theory, the additional scalar 
field called quintessence can explain the current acceleration of the Universe^^^. 
(II) Inflation can be explained in the context of scalar-tensor theory. This inflation 
model is called the hyper- extended inflation^'^® where the additional scalar field acts 
as the inflaton. (HI) Scalar-tensor theory also appears in the particle-physics. Su- 
perstring theory^'^^^ predicts that the dimension of the spacctime is 10 or 11, and 
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in the low energy effective theory (supergravity theory), scalar fields called dilatons 
appear. In this theory, the gravitational theory is of the scalar-tensor type. Also, in 
Kaluza-Klein theory where GR is extended to 5 dimension with the extra dimension 
being compactified, it reduces to the scalar-tensor type with 355 component acting 
as a scalar field 

As for GWs from compact binaries, there exists scalar dipole radiation which 
modifies the binary evolution from GR. After introducing the action in Sec. 6.1, 
we derive the correction to GW phase in Sec. 6.2. Then, in the final section, we 
describe the current constraints from solar system experiment and binary pulsar 
observations. 



6.1. The Action 

In general, the action for the scalar-tensor theory in the Jordan frame ^ is given as 

1 ..1 



^ " 16^ i d -^V^ [4>R - uj{c^)^g^''d^(j)d,ct) + 16^L„,attcr(*) j ■ (140) 

Here, is the scalar field that carries the scalar gravitational degrees of freedom 
and 4" is the matter field. The first term is the Einstein-Hilbert term non-minimally 
coupled to the scalar field. The second term is the kinetic term of the scalar field. 
The coefficient w((/)) is the parameter that characterizes this theory. Especially, when 
a;((/)) = wbd = const., this theory is called Brans-Dicke theory and wbd is called 
Brans-Dicke parameter. (The inverse of) this parameter expresses the strength of 
the coupling between the scalar and the matter fields (c.f. Eq. (147)). If we take the 
limit wbd ^ oo, the theory reduces to GR. 

The third term in Eq. (140) is the Lagrangian for the matter field. At the micro- 
scopic level, this term does not depend on <f>. However, if one describes the orbital 
dynamics of a system of strongly self-gravitating objects by means of an effective 
point-particle theory, in which the internal structure of the objects is integrated out, 
then ^matter consists of a sum of point-particle Lagrangians which do depend on 
4>, on account of strong equivalence principle violation. For example, see Eq. (145) 
below. 

In this theory, the gravitational constant is not a constant but a function of (p. 
It is given as^^^ ^ 

GBD(0)^i±|^^. (141) 
If we assume that the theory reduces to GR at spatial infinity, we have 

4 + 2WBD 1 



Gbd{x oo) = 



3 + 2wbd 00 

G = l, (142) 



•^Somc literature uses a different conformal frame (the Einstein frame)^^^' 

'This can be derived by perturbing the scalar and the metric fields around the flat background to 
obtain the gravitational potential. 
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where 0o is the value of the scalar field at spatial infinity. Therefore, we obtain 

where 7 is a parameterized post- Newtonian (PPN) parameter^^'' defined as 

7^^. (144) 

In GR, 7 takes the value 7 = 1. 



6.2. Correction to GW Phase of the Binary 

In this section, we derive the correction to the gravitational waves from compact 
binaries in Brans-Dicke theory, assuming that wbd ^ 1- First, the action for the 
binary system in this theory is given as^^^^ s 

^ " i / '^'^^ {f^ ~ ^'^■"'^'") " ^. / ^^^^^ 

where Ta is the proper time measured along the trajectory of a-th body. The field 
equations are given as 

Gfj.!^ ^^^"^^"^ + - ^^M'^'/'^Ai'/''^^ + ^('?^;m>' ~ 9f^iyO(f>), (146) 

where G^^ = R^^i, — ^Rg^i, is the Einstein tensor and Tfj^n is the matter energy- 
momentum tensor. The scalar wave equation can be obtained as 



dT 

r-20— 



(147) 



3 -I- 2wbd 

In the point-particle limit, T^ii^ and dT/dcj) can be expressed as 

T^" = ^ E rna{c^)'^5\x - x^), (148) 



^ a=1.2 



ST_ 1 ^5^1,3,, ,,,,, 



Next, let us consider the perturbation of g^u and (j) around spatial infinity: ^^"^ 
9tiy = ilfiu + h^„, (/) = (t)Q + (t>. (150) 
We define a new tensor field 

Xf^i' = hf,^ - l;hr]f,^ - -^rif,^, (151) 
Z (po 



sSee also Refs.^'*"' ^''^ where the authors put the point-particle action into the modern context of 
Effective Field Theory, and interpreted rua {<j>) as the leading-order term in the derivative expansion 
of the point-particle Lagrangian of body a. 
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and impose the Lorentz gauge condition 

By perturbing the field equations (Eqs. (146) and (147)), we get 



□0 = -leTrS", 



(152) 

(153) 
(154) 



where the box operator represents the d'Alenibertian operator in Minkowski back- 
ground and S is defined as 



S 



1 



T-200 



dT 



(155) 



6 + 4wbd 

where the index denotes the value at spatial infinity. 

At a distance r that is sufficiently larger than the size of the binary system, we 
can perform the multipolar expansion to yield^^"* 



X 



d 

ml \ dt 



r ^ m! V 

m=0 



S(t-r, x'){n-x'Y^d^x' 



(156) 



(157) 



First, we consider the correction to the tensor part. From the gauge condition and 
the field equations, we have T^^ — 0, and hence at the leading order. 



r <pQ 
1 + 7 cf^ 



r dt'^ 

where we have used Eq. (143). At the lowest order, we have T^'^ « p so that 

1+7 



dt^ 



Ei j 



(158) 



(159) 



a=l,2 



The term proportional to 7 represents the correction compared to GR. However, 
we now show that this OPN correction is sub-leading compared to the one from the 
scalar field. 

Next, let us consider the contribution from the scalar field. First, we define the 
sensitivity (which is related to the scalar charge) as 

9(ln ma) 



Then, S can be rewritten as 



1 



6 + 4wbd 



9(lnG) 



(160) 



(161) 
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Wc substitute this into the scalar wave equation and solve for (f). At spatial infinity, 
we can assume that 0o is conserved. This means that the mass is conserved, and 
hence the leading order scalar monopole radiation vanishes. Therefore the leading 
radiation appears at dipole order. When we move to the center of mass frame with 



(162) 



the scalar perturbation (f) can be expressed as 



-Sn-v (163) 



3 + 2a;BD r 

with S = S2 — si. This S in Eq. (163) comes from "-2sa" part in Eq. (161). The 
first term in Eq. (161) is sourced by the mass of each component, but because of 
the linear momentum conservation, there is no dipole radiation from this term. 

Next, we calculate the correction to the radiated energy In order to perform 
this, we need to expand perturbation up to the second order. Then, Eq. (153) 
becomes 

Ux^" = -167rT^'^, (164) 

where 

r^'^ = -^T^'^ + t^"" . (165) 

00 

Here, t'"^ is a second order tensor composed of and cj) corresponds to the 
energy-momentum tensor of GWs. (The explicit form is given in Appendix A. 3.) By 
taking the TT gauge and performing integration by parts, we obtain the expression 
for the GW energy flux as 

dt " 



327r 



dAt" 

00 / / ixli'xli' + (4u;bd + 6)^^]dn ) (166) 



with (fi = 0/00- By substituting Eqs. (159) and (163) into the above equation and 
keeping only the leading PN order for the correction, we obtain 



— ' -V 1 + —o ujv ) (167) 



dt 5 



1 + —S 



2 - -2 
UJV 



with a) = Wgj!). The second term corresponds to the correction due to the scalar 
dipole radiation. The monopolar radiation (not included in the above equation) 
appears at OPN order relative to GR,^'^^ which is subdominant compared to the 
dipolar radiation. 
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Finally, we derive the correction to the gravitational waveform phase"'^^^. By 
using Eq. (167), the time evolution of frequency can be calculated as 

I _tl _ I dVl da dE 
TT TT da dE dt 



3 Af//^ 2a2 



Egw 



96 1 ^A/A"/' / , 5 _i 



By performing the integration once, we derive the correction to t(/) as 

t{f) = to - 5X(87rX/)-8/3 (^1 _ ^S^Ojx-^^ , (169) 
and the gravitational waveform phase in the Fourier domain is given as ^ 

*(/) = 27r/to - 00 - 7r/4 + *opn(/) (^1 - ^S^ux-^^ . (170) 

Higher PN terms in GR that appear in Eq. (96) can be added to this equation 
accordingly. 

From the corrections above, we see that the effect of the scalar dipole radiation 
is greater for binaries with larger S. The sensitivity Sj, roughly corresponds to 
the self-gravitating energy of the body with svvd ^ 10~'^, sns ^ 0.2 and sbh = 
0.5 for WD, NS and BH, respectively i^^. Therefore, binaries whose components 
consist of different types of compact objects give larger BD effect. In this theory, 
the (rescaled) scalar charge is given as q = 1 — 2sa- Hence, BHs do not have scalar 
charges and this is a consequence of the no hair theorem^*'^"^*® that also holds 
in BD theory^^^ However, the proof of the no- hair theorem relies strongly on 
the assumption of stationarity, and there are astrophysically interesting situations 
where this assumption is violated and BHs grow scalar hair'^^-'^'^^^. 



6.3. Current Constraints 
6.3.1. Solar System Experiment 

In this subsection, we explain the current constraint on BD theory obtained from 
the Shapiro time delay measurement by Cassini satellite^^°. First, the metric for 

the solar system can be written in PN gauge as^^*"' 

ds^ = -(1 - 2$)dt2 + (1 - 2-f<i>)6,jdx'dx^ + ©(e^/^) , (171) 



''See Ref.^*^ for higher PN corrections at linear order in r} obtained by applying Teukolsky formal- 
ism. Also, very recently, Yunes et al}^^ solved the Teukolsky equations numerically and obtained 
a fitted PN formula. 

'Sec a related work by Sotiriou and Faraoni^^*^ where they extend Hawking's result to more general 
class of scalar-tensor theories. 
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where the dimensionless parameter e is defined as e = Mq/v with r denoting the 
typical length scale and $ = —AIq/t + 0{e^). 

The most stringent constraint on 7 has been obtained by Cassini spacecraft 
while on its way to Saturn. The light from the satellite that passes close to the 
sun gets deflected and there arises a time delay called Shapiro time delay. (See e.g. 
Ref.^^'^ and references therein.) We take this path of the light to be a straight line 
parameterized by a; = (x, b, 0) where b is the impact parameter and 
for the satellite and the earth, respectively. For a round trip, the additional time 
delay can be expressed as 

At ^-2 (i+j)^^^^-^dx, 

= (l+7Kln^, (172) 

where rg = 2Mq is the Schwarzschild radius of the sun, a® = lAU and = 8.43AU 
is the distance from the satellite to the sun. Here, we have assumed that ^ 
b and similar condition for the satellite. Cassini measured the frequency shift 

_ dAt _ dAt db 

~ ~dr ~ ~db"dt ' 
4r„ db 

and obtained the result that (1 + 7)/2 must be within 0.0012% of unity. This gives 
the constraint J 

CJBD > 4 X 10'' . (174) 

This constraint is expected to be improved by 3 orders of magnitude using AS- 
TRODI61. 



6.3.2. Binary Pulsar Test 

There are also stringent constraints obtained by the binary pulsar tests. As we 
explained in Sec. 6.2, there exists a scalar dipole radiation in BD theory, and hence 
we can place constraints on the theory from the observations of the orbital decay 
rate of the binary pulsar systems. All accurately-timed binary pulsars place a bound 
on cjbd as shown in e.g. Fig. 7 of Ref.^^^. (See also Ref.-'^^^.) Among them, the binary 
pulsars that consist of a NS and a WD perform better in constraining the theory 
compared to the ones consist of 2 NSs, for the reason explained at the end of Sec. 6.2. 
The WD-pulsar systems PSR J1141-6545 and PSR J1738+0333 place the 
bounds as 

r2xl04 (J1141 - 6545), 
WBD > < . ^ ^ (175) 

[2.5x104 (J1738 + 0333), 



■■See related works by Perivolaropoulos'^^* and Alsing et al.^^^ for the current constraints on 
massive BD theories. 
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respectively. These are slightly weaker than the one obtained from the solar system 
experiment, but still they are meaningful in the sense that they constrain the dipole 
radiation of the binary in BD theory. 

7. Massive Gravity Theories 

Massive gravity theories are simple extension of GR where we add finite mass to 
graviton. There are many different types of massive gravity theories (for recent 
reviews, see e.g. Refs."'^^'^'^^^). Originally, Fierz and Pauli"'^''^ proposed a Lorentz- 
invariant massive gravity by simply adding a quadratic graviton mass term to the 
Einstein-Hilbert action. However, it was found that the linearized Fierz-Pauli the- 
ory does not reduce to the linearized GR in the massless limit (van Dam-Veltman- 
Zakharov (vDVZ) discontinuity^^'"' ^^^). If this discontinuity is valid in the solar 
system, Fierz-Pauli theory is ruled out by the solar system experiments. How- 
ever, Vainshtein proposed that the effect of nonlinearity cannot be neglected^^^ and 
unlike GR, linear approximation breaks down already at a distance much larger 
than the Schwarzschild radius (Vainshtein radius). Later, Nicolis and Rattazzi^^^ 
showed that this Vainshtein mechanism indeed works for the DGP braneworld 
model^^''. Rubakov^^^ and Dubovsky^^^ proposed Lorentz- violating massive gravity 
theories which evade pathologies related to the vDVZ discontinuity. Chamseddine 
and Mukhanov came up with a massive gravity model that is analogous to the Higgs 
mcchanism"'^^'^. They introduced 4 scalars with global Lorentz symmetry. When it 
is broken, graviton absorbs scalar degrees of freedom and acquires a finite mass. 
Although Vainshtein mechanism seems to work in this theory^^^, it is realized that 
the appearance of ghosts cannot be evaded^^^. 

Recently, de Rahm et al}^^''^^'^ proposed a novel massive gravity (the so-called 
non-linear massive gravity) under a flat reference metric, which is a non-linear gen- 
eralization of Fierz-Pauli theory. In Ref.^*^, it has been shown in the unitary gauge 
that the Hamiltonian constraint which is required to kill BD ghost exists up to 
fourth order in non-linearities for a certain range of free parameters. Later, this 
theory is generalized to generic reference metric by Hassan and Rosen^^'^. Then, 
the existence of the Hamiltonian constraint has been shown to exist to all orders 
in the unitary gauge with a reference metric taken as (i) flat^^*, (ii) generic but 
non-dynamical, and (iii) generic and dynamical (non-linear bimetric gravity)^^^. 
Moreover, the secondary constraint is also shown to exist^^^ for both non-linear 
massive and bimetric gravities, concluding that BD ghost is evaded up to all or- 
ders in these theories. See Refs.^^®'^^^ for the proofs of existence of the secondary 
constraint in generic gauge. 

Current constraint from solar system experiment puts a model-independent con- 
straint as Xg > 2.8 X lO-'^^cm As for GWs, the propagation speed is modified 
from the speed of light when graviton acquires a finite mass. In the next subsection, 
we derive the correction to the GW phase due to this deviation in the GW prop- 
agation speed from c. Then, in Sec. 7.2, we explain the current constraints on the 
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Fig. 13. Observer detects a GW signal at t = fa that is emitted from X = Xe at fe. Then, he 
detects another signal at t = that is emitted at t = t'^. 

mass of graviton both from the solar system experiment and binary pulsar tests. 
7.1. Correction to GW Phase of the Binary 

When the mass of graviton is non- vanishing, the propagation speed of GW becomes 
less than the speed of light, which shifts GW phase from GR. In this section, we 
derive this correction to GW phase based on Will^®^. 

Let us first assume that the graviton propagates in a homogeneous and isotropic 
Friedmann-Lemaitre-Robertson- Walker spacetime. Its metric is given as 

ds^ = ~dt^ + a'^{t) [dx"^ + J:^{x){d9^ + sin^ 9d(P^)] , (176) 

where a{t) is the scale factor of the background universe and S(x) equals sin^, X 
or sinh x for spatially closed, flat or open universe, respectively. When a graviton 
is radially propagating from x — Xe the detector % = 0, we can show from the 
geodesic equation that = const., where is a 4-momentum vector^'^^. The norm 
of p'^ is given as 

P^'p'^g^u = -E' + ^^ -ml , (177) 
where E = p^ . From this equation, we obtain 



-1/2 
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with 

pl^a\U){El~ml), (179) 

where is the emitted energy of the graviton. Let us assume that 3> mg and 
expanding Eq. (178) to first order in {mg/Eef' ■ After taking integration, we obtain 

dt 1 ml /■*" 

-TT" 9 2(, \j,2 / "W'^^' (180) 
a{a) 2a^{te)E^Jt^ 

where ta is the arrival time of graviton at % = 0. Now, we consider two gravitons 
emitted at tg and t'^ with energies Ee and E'^, and arrive at ta and f^, respectively. 
(See Fig. 13.) Assuming that Ate = te ~ t'^ ^ a/a, we obtain 

D 



Ata ^{I + Z) 



(181) 



a 

where Ata ^ ta — t'^ and fe is the emitted frequency which can be expressed as 
/e = Ee/{mgXg). The new distance parameter D is defined as 

D = ii±il /*° a(Odt 

with tto = a(ia) representing the present scale factor. 

Under the stationary phase approximation, the phase of GW in Fourier domain 
is written as 

*(/) = 2tt [\t - to)df + 2TTfto - 00 - 7 , (183) 

where / is the observed frequency and /o is the one at coalescence. By substituting 
Eq. (181) into the above equation and using /e = (1 + z)f, we have 

*(/) = 2tt / (te - teo)dfe - TT2 + 2^/^o - 00 - 7 , (184) 

where t^o and /eo are the coalescing time and frequency measured in the source 
frame and 

= 2(1 + .)Xlfr '^° = ^°-(i + .)A^/o- ^'''^ 

By performing integration and re-expressing Eq. (184) in /, we obtain 

^if) = 27r/to - 00 - J + *opn(/) - ^I3g7j'/'x^ , (186) 

where 

(When we perform Fisher analysis later, we replace to and 0o with ^o and 0Oj 
respectively.) 
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7.2. Current Constraints 

In this section, following Ref.^^^, we explain the current static and dynamical con- 
straints on the mass of graviton. We first explain the current static constraints 
mainly from solar system experiment When the mass of graviton is finite, 
the gravitational potential is given as a Yukawa type: 

V{r) = -^c^p(^-^y (188) 

Under the assumption r ^ Xg, the gravitational acceleration is estimated by 
taking the derivative of the above equation with respect to r as 

(189) 

where n* is the unit radial vector and 

<;„._(,,i_)„p(_i_) 

From the observation of the Kepler's third law (a(27r/P)^/^ = G{a)^^^ with P 
denoting the orbital period) for Mars, the bound on Ag has been obtained as 

Ag > 2.8 X 10"cm. (191) 

Other static constraints include the evidences for the bound clusters and tidal 
interactions between galaxies^™, and the bound from weak lensing^''^. Although 
these bounds are considerably stronger than the solar system bound, there remain 
large ambiguities in these bounds due to the uncertainties in the amount and dy- 
namics of dark matter. Recently, Sjors and Mortsell^''^ calculated the constraint on 
Ag in non-linear massive gravity^^^' -'^^^ using galactic lensing and velocity dispersion 
data. Under the assumption of static, spherically symmetric spacetime, the bound 
is given by Ag > lO^^cm. We stress here that this is a model-dependent constraint 
while the one from the solar system experiment does not rely on the specific choice 
of the massive gravity theories. 

Other than the static constraints, there is a dynamical one obtained from a 
binary pulsar test^''^. For a specific type (Fierz-Pauli type) of the massive gravity 
theory, they estimated how radiated GW luminosity is modified in a binary pulsar, 
and obtained a constraint from PSR B1913+16 and PSR B1534+12 as 

Ag > 1.6 X lO^'^cm, (192) 

which is 2 orders of magnitude weaker than the solar system experiment. These 
current constraints are summarized in Table 2. 
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Table 2. The current constraints on X„. 



Current Bound 


Ag(cm) 


Solar System^*'' (model-independent) 

Clusters-^^" (model-independent) 
Weak Lensing-^''^ (model-independent) 
Galaxies^''^ (Nonlinear Massive Gravity) 
Binary Pulsars^'^^ (Fierz-Pauli type) 


2.8 X 10^' 
6.2 X Vp-^ho 
1.8 X lO^'^ 

1026 
1.6 X 10^^ 



10^ 



10' 



eq.-mass (DPF) 
eq.-mass (adv. DPF) 
IMRI (adv. DPF) 
SNRthr. 




10° 10^ 10^ 10^ lO"^ 10^ 10^ 10'' 
time to coalescence (sec) 



Fig. 14. Accumulated sky-averaged SNRs with DPF (red thick solid) and with adv. DPF (green 
thin solid) against the time to coalesce for a 10^ Mq equal-mass IMBH binary with x = 0.2 in NGC 
6752. We also show the one for a (10 -|- 2 X W^)Mq binary in NGC 6752 with adv. DPF (blue 
dotted). The horizontal dotted-dashed line corresponds to the detection threshold of pthr = 5. 
(This figure is taken from-*^®*.) 



8. Observing GWs from IMBHs with DPF 

Although the event rate is not so promising, in this seetion, we consider how DPF 
might help ground-based detectors as a complementary observation and how accu- 
rately one can measure binary parameters if the signal is detected. In this section, 
we do not consider deviations away from GR. 



8.1. Possible Joint Search with Ground-Based Detectors 

From Fig. 4. one sees that the GW signals from an equal-mass IMBH binary in NGC 
6752 can be detected by both DPF and ground-based interferometers, and hence 
the joint searches of these signals with these detectors may be possible. To clarify 
this possibility, we show the accumulated SNRs of DPF and adv. DPF against time 
to coalesce for an equal-mass IMBH binary and an IMRI in Fig. 14. Let us focus 
on the former case, where the accumulated SNR reaches p = 5attw4x lO^s. 
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(The definition of SNR is given in Eq. (195).) This means, in principle, DPF will 
be able to give an alert to ground-based detectors about 7 mins before coalescence 
to make sure that they are operating at the time of merger. By combining the DPF 
and the ground-based data, we would be able to observe the whole history of the 
late inspiral, merger and ringdown. (Sec Rcf.^''^ for a related work on possible joint 
observations of IMBH binaries with LISA, ET and adv. LIGO.) DPF data may help 
in confirming GW signals for the ground-based detectors when only the ringdown 
signals have been detected. Notice that the times the accumulated SNR reaches 
p = 5 are almost identical for DPF and adv. DPF. This shows that it is more 
important to improve the sensitivity at lower frequency to claim earlier detection. 
For example, if we improve DPF sensitivity by a factor of 2, the accumulated SNR 
would reach p = 5 about 1 hour before coalescence. 



8.2. Fisher Analysis 

We use the matched filtering analysis to estimate the determination errors of the 
binary parameters 6 26,275 assume that the detector noise is stationary and 
Gaussian. The noise takes the Gaussian probability distribution given by 



p{no) oc cxp 
where we have defined the inner product as 



{A\B) = 4Re / df 



A*{f)B{f) 



Snif) ■ 

The signal to noise ratio (SNR) for a given GW signal h is 



p[h] ^ ^/{h\h). 



(193) 



(194) 



(195) 



The detected signal s{t) is the sum of the gravitational wave signal h{t; Of) and 
the noise n{t), with dt representing the true binary parameters. Then, Eq. (193) 
can be rewritten as 



p(0t|s) cxp(o)(0t)exp 



(ht\s)~^{ht\ht) 



(196) 



where p^^\Ot) denotes the distribution of prior information and ht = h{6t). We 
determine the binary parameters as 6 that maximizes the probability distribution 
p{6t\s). 6 satisfies the following equation. 



(197) 



with di = . We can express 9^ as 0* = ^* + with representing the error in 
the determination of 6'\ Next, we expand Eq. (196) in powers of up to quadratic 
order and find 



p{e\s) cxp(°)(0)cxp 



1, 



) 'j 



(198) 
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estimation error 



6' 



Fig. 15. Gaussian probability distribution of binary parameters 6' 



where = {didjh\h — s) + {dih\djh) is called the Fisher matrix (see Fig. 15). Since 
h — s = —n, we can neglect the first term of Tij in the limit of large SNR and 



r„- = {d,h\djh). 



(199) 



When we perform sky-avcragcd analysis under the assumption that N-mt effective 
interferometer arc placed on one site, the Fisher matrix can be written as 



r,, = Ni^tidrh\djh). 



(200) 



We set the prior information that the dimensionless spin has to be less than 1. 
This can be achieved by setting^^^ 



3(°)(0) ocexp 



2 I 9.4 



1 / cr \2 



2 V2.5 



(201) 



The root-mean-square of A9^ can be calculated by taking the square root of the 
diagonal elements of the covariance matrix E'-^ , which is the inverse of the Fisher 
matrix: 



(202) 



with Tij defined by 



p(°)(0) exp 



2 ' 



cxp 



-^f,jA9'Ae^ 



(203) 



We take the integration range of the inner product in Eq. (194) as (/in, /fin) 
where 



/in = max{/iow, /tobs}, /fin = min{/high, /isco}- 



(204) 
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/tobs is the frequency at tobs before the binary reaches ISCO and it is given by 

-5/8 / , \ -3/8' 



5.5 X 10" 



M 
lOMp 



( ^obs 



Hz. 



(205) 



Here, we only take the leading contribution from GR into account. The frequency 
at ISCO is given by 



(206) 



We recall that /low and /high correspond to the low and high cutoff frequencies of 
the detector, respectively. 

8.3. Binary Parameter Estimation 

In this subsection, we explain how accurately DPF can determine binary parameters 
by using the inspiral-merger-ringdown hybrid waveform. 




10' 10-' 




Fig. 16. (a) DPF Parameter determination accuracies of In (red thin solid), Inr; (green thick 
solid), X (blue thin dotted) and XnD^ (magenta thick dotted) against the total mass Mt for equal- 
mass binaries of x = 0.2 at = 5kpc, obtained from pattern-averaged analysis. The vertical 
dotted-dashed line represents a possible IMBH binary in lo Centauri. (b) The SNRs for equal-mass 
binaries of x = 0.2 at = 5kpc against the total mass Mt with DPF (red thick solid) and with 
adv. DPF (green thin solid). The horizontal dashed line corresponds to the detection threshold 
Pthr = 5. (This figure is taken from Ref. ^^''.) 



8.3.1. Pattern-Averaged Analysis 

First, we show the results for the Pattern-averaged analysis, where we take the 
sky-average of the waveform. There are 6 parameters in total: 

0* = (lnMt,ln77,x,to,0o,^L). (207) 
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Fig. 17. (a) Histograms showing the parameter determination accuracies of In Mt (red thin sohd), 
In Tj (green thick soUd) and x (blue dotted) using DPF obtained from Monte Carlo simulations. We 
assumed 10'' Mq equal-mass IMBH binaries in uj Centauri with x = 0-2 and random orientations 
for the orbital angular momenta, (b) Histograms showing the parameter determination accuracies 
of InDi (magenta dotted-dashcd) and Qg (blue solid) using DPF. (This figure is taken from 
Ref.W4.) 



We take the fiducial values of mi/m2 = 1, x = 0.2, = 0o = 0, and Dl = 5kpc. 
The panel (a) of Fig. 16 shows the DPF measurement errors of InM* (red thin 
solid) , In -q (green thick solid) , x (blue thin dotted) and In (magenta thick dotted) 
against the total mass Mt- The dottcd-dashed vertical line corresponds to a possible 
equal-mass IMBH binary in uj Centauri. One can see that DPF has an ability to 
measure binary parameters within several % accuracies. The corresponding SNRs 
for DPF and adv. DPF are shown in the panel (b). The parameter estimation for 
adv. DPF can be roughly estimated by linearly scaling the DPF results with the 
SNR ratio. For example, the ones for an equal-mass IMBH binary in w Centauri 
improves roughly by a factor of 7. Since the SNRs are not so high, the Fisher analysis 
can only give rough estimates^^^. 



8.3.2. Monte-Carlo Simulations 

Next, we perform Monte-Carlo simulations, where we randomly distribute the di- 
rections and orientations of the sources. This time, we have 10 parameters in total: 

= {\nMtMTj,x,to,(f>Q,DL,Os,^s,ei^,^L), (208) 

We set {mi, 7712) = (lO'*, 1O^)M0, x — 0.2 and Dl = 5kpc. We use the galactic 
latitude and the longitude of w Centauri, which can be converted into {9s, 4's) using 
the formulas given in Appendix A of Ref.^^^. We randomly generate 10"* sets of 
COS0L in the range [-1,1] and (pL in the range [0, 2tt]. We also distribute (fEo and 
t^Do (the detector position at t = 0) in the range [0,27r]. For each set, we calculate the 
Fisher matrix. 10"* sets of results allow us to construct the probability distribution 
of each binary parameter. 
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In the panel (a) of Fig. 17, we show the probabihty distributions of In Afi, Inry and 
X- One sees that these parameters can be determined within several % accuracies. 
In the panel (b), we show the results for InD^ and fts where the latter is defined 
as 



Ails ^ 2n\ sin0-s|y'E,,,,S]^3^3 - S]|^^^. (209) 

Unfortunately, these parameters cannot be determined with high accuracy. The 
angular resolution can be roughly estimated as 

/1A\' /30\V A \ V6.4X lO^cmA' , , 

Here, A is the wavelength of GWs and D is the effective size of the detector which 
corresponds the diameter of Earth for DPF observation. This estimate agrees with 
the Monte-Carlo results. Dl is determined from the amplitudes of GWs, but since 
the angular resolution is of 0(1), InD^ is only determined to the same accuracy. 
This shows that distance and angular parameters are poorly determined by DPF. 
However, because the correlation between the amplitude and the phase parameters 
are very weak, DPF would be able to determine the mass and spin parameters with 
good accuracy. Since these signals cannot be detected by ground-based detectors, 
DPF would perform unique observations. 



9. Testing Alternative Theories of Gravity v^rith Gravitational 
Waves 

One of the most interesting scientific goals of GW observations is to test alternative 
theories of gravity in the strong-field regime. In this section, we review the proposed 
constraints on Brans-Dicke and massive gravity theories by observing gravitational 
waves from compact binaries. 



9.1. Brans-Dicke 

9.1.1. Proposed Constraints with ground-based detectors and LISA 

Constraining Brans-Dicke theory using gravitational waves from compact binaries 
was first discussed by Eardley-'^^^. In this theory, the non- vanishing scalar dipole 
radiation exists^^^' -"^^^j which modifies the binary orbital evolution from the one 
in GR. The change in the orbital evolution due to this dipole radiation modifies 
the phasing of the gravitational waveform. Will^^^ carried out the matched filter- 
ing analysis and calculated how accurately one can measure the binary parameters 
including lubb with adv. LIGO. Damour and Esposito-Farese^^^ extended Will's 
work to generic scalar-tensor theories. Scharre and Will-^^'' and Will and Yunes^^® 
followed Will's analysis using LISA. Will and Yunes improved the previous works 
by investigating how the constraints depend on LISA position noise, acceleration 
noise and arm lengths. These works consider non-spinning binaries under pattern- 
averaged (sky-averaged) analysis. However, these calculations do not include binary 
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spins and also they used the pattern-averaged waveforms. Berti et al}^'^ calculated 
the determination accuracies using LISA, taking the effect of the spin-orbit coupling 
into account for spin-aligned binaries. They performed both the pattern-averaged 
and the Monte Carlo simulations. In the latter simulations, they randomly dis- 
tribute the directions and the orientations of IC^ sources over the sky, calculate 
the constraint on wbd from each binary and take the average at the end. These re- 
sults are summarized in Tables 3 and 4. There are also some works on constraining 
massive BD theoryi^'i'i'^^ ^nd generic theories that produce dipolar gravitational 
radiation^^"* with GW interferometers. Very recently, merger simulations of BH^^'^ 
and NS^^* binaries have been performed. 

Tabic 3. Summary of the works on probing BD theory using GWs from compact binaries. The 
second column denotes up to which PN order they take. The third, fourth and fifth columns represent 
whether they take the spin-orbit coupling /3, precession or eccentricity into account. The sixth column 
describes whether they assume multiple-source detections. The seventh column shows what kind of 
analyses they have performed (either the pattern-averaged (PA) or Monte Carlo (MC).) Finally, the 
last column denotes whether generic scalar- tensor theories are considered. 



Reference 


PN 


/3 


prec. 


ecc. 


multi. 


analy. 


generic ST 


Will (1994)^'"'' 


1.5 


X 


X 


X 


X 


PA 


X 


Damour & Esposito-Farese (1998)^^^ 


1.5 


X 


X 


X 


X 


PA 


o 


Scharre & Will (2002)167 


1.5 


X 


X 


X 


X 


PA 


X 


Will & Yunes (2004) ^i^^ 


1.5 


X 


X 


X 


X 


PA 


X 


Berti et al. (2005) ^'^^ 


2 


o 


X 


X 


X 


MC 


X 


Yagi & Tanaka (2010)"° 


2 


o 


o 


o 


X 


MC 


X 


Yagi & Tanaka (2011)"^ 


2 


o 


o 


o 


o 


MC 


X 



Table 4. Summary of the proposed constraints on wbd using GWs. The numbers in the brackets 
denote the total masses of the binaries in the unit of Mq . 



Reference 


adv. LIGO 


LISA 


DECIGO/BBO 


Will (1994)168 
Damour & Esposito-Farese (1998)1^'^ 
Scharre & Will (2002)1^7 
Will & Yunes (2004)168 

Berti et al. (2005)169 
Yagi & Tanaka (2010)"'' 
Yagi & Tanaka (2011)"^ 


2 X 10^ (3.7) 
2.2 X 10^ (11.4) 


2 X 10^ (10^) 
2 X 10^5 (10^) 

lO* (10^) 
7 X 10^ (10^) 


4 X 10® (11.4) 



In Ref.^^°, we further improved these analyses by (i) taking the spin-spin cou- 
pling effect into account, (ii) considering slightly eccentric binaries, and (iii) looking 
at precessing binaries^^''. We first performed the pattern-averaged analysis where 
we have 9 parameters in total: 

9' = (ln7W,ln7/,/3,a,to,0o,i?L,/e,w). (211) 

We set the fiducial values oil3 = a = iQ = (f)~iI)~Q and set p = 10. For the 
initial eccentricity, we assume — 0.01 at 1 yr before coalescence. We consider 
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BH/NS binaries and set S = 0.3. Wc fix mNs = 1.4M0 and consider 4 different BH 
masses, tobh = 400, 1000, 5000 and lO"' Mq. For each binary, we show the results 
for circular and eccentric binaries in Table 5. (a is included in both cases.) For 
comparison, we also show the results in Ref.^^^ where the authors did not include 
cr nor le into binary parameters. 

From this table, one can see that the determination accuracies deteriorate as 
we increase the number of parameters. This is because the parameters are strongly 
correlated and adding parameters dilutes the binary information in the detected 
GWs. Inclusion of both a and /g into parameters increases the determination errors 
by roughly 1 order of magnitude. In particular, the inclusion of has larger effect on 
determining wbd than the inclusion of a. The reason can be understood as follows. 
In the phase ^'(/), the term containing wbd has a relative frequency dependence of 
/~^/^ compared to the leading Newtonian term, while the ones containing a and le 
have the relative frequency dependences of f^^^ and f~^^^^, respectively. The first 
(BD correction) and the third (eccentricity) terms both have negative power-law 
indices which makes the correlation between these 2 parameter stronger than the 
one between wbd and a. 

The table also shows that the constraint becomes more stringent as the BH mass 
decreases. This is because the orbital velocity of the binaries become smaller, which 
makes relative "-1PN" dipole correction larger. In Ref.^™, we also show that the 
measurement accuracy of wbd would be similar to the one without including 1^ if 
we impose the prior information of /g > 0. 

Table 5. The results of error estimation in Brans-Dicke theory using LISA^™. These are calculated 
with pattern-averaged analysis for the spin-aligned BH/NS binaries with various BH masses under 
p = 10. We fix NS masses to 1.4Mq. For each binary, the first line represents the results obtained 
by Berti et al.^^^ which does not include a nor le into parameters. The second line represents the 
results including cr, while the third line shows the ones including both cr and le- We used only one 
interferometer for the analyses. (This table is taken from Ref. ^^''.) 



cr 7e 




AlnA^(%) 


Alnr; 


A/3 


Ato{s) 


Mo 


Ao- 


A/e(10--^") 


4OOAf0 


















X X 


39190 


0.00657 


0.0250 


0.0508 


7.95 


76.7 






O X 


24886 


0.0130 


0.0819 


0.202 


13.8 


552 


2.39 




o o 


4583 


0.0396 


0.142 


0.280 


16.7 


552 


2.49 


1.09 


lOOOA'/© 

X X 


21257 


0.00764 


0.0186 


0.0557 


7.99 


58.4 






O X 


8210 


0.0265 


0.110 


0.0692 


23.5 


919 


1.96 




o o 


1881 


0.0692 


0.193 


0.261 


23.6 


1059 


2.41 


6.34 


5OOOM0 

X X 


6486 


0.0114 


0.0133 


0.0550 


8.79 


23.4 






O X 


1933 


0.0503 


0.0936 


0.221 


37.9 


1108 


0.595 




o o 


281 


0.224 


0.302 


0.916 


62.9 


2438 


1.30 


173 


IOOOOMq 

X X 


3076 


0.0178 


0.0161 


0.0706 


13.6 


15.5 






O X 


862 


0.0827 


0.114 


0.350 


82.9 


1763 


0.474 




o o 


113 


0.412 


0.418 


1.51 


160 


4454 


1.11 


797 
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Next, we show the results for the Monte Carlo simulations. We have 13 param- 
eters in total: 

9' = (lnAl,ln77,/3,a,<o,0o,i?L,/e,^s,0S,^L,0L,w). (212) 

We fix the masses of the binary constituents as tons = 1-4A/0 and tobh ~ 10'* M©, 
and consider a binary of p = \/200. (This roughly corresponds to p = 10 for a single 
interferometer). We randomly distribute eos^s and cos^l hi the range [-1,1], and 
and in the range [0,27r], generating 10"* sets of spin- aligned binaries. We calculate 
the measurement accuracies of w for each binary using Fisher analysis and find the 
probability distribution of the future lower bound on wbd shown as the (light blue) 
thin dotted histogram in Fig. 18. For comparison, we also show the solar system^^° 
and the binary pulsar^^^ constraints as the dashed and the dotted-dashed lines, 
respectively. This histogram shows that LISA can only place weaker bounds than 
the current ones. Notice that the histogram has "tails" on both weaker and stronger 
bounds. The former corresponds to binaries with L ■ N ~ Q. The latter is due to 
the Dopplcr phase 

(pz,(t) = 27r4^ sin 0s cos[,^(t) - ,^s] , (213) 

where 

fc=^ = 2.00mHz , (214) 
K 

is the critical frequency with R = lAU. This Doppler phase makes the degeneracies 
between some parameters stronger, but for some specific directions, the effect of the 
Doppler phase vanish and the constraint becomes stronger. For example, the deriva- 
tive oi ^pD with respect to is proportional to cos^Si a-nd hence the correlation 
between wbd and due to the Doppler phase is disentangled at = 7r/2. 

Next, we show the results for the precessing binaries. There are 15 parameters 
in total: 13 parameters for the spin-aligned case plus n and Uc- For Monte Carlo 
simulations, we set the dimensionless spin parameter of a NS and a BH to and 0.5, 
respectively. We randomly distribute the fiducial values oi k = L ■ S ui the range 
[-1,1] and ttc in the range [0,27r]. 

The results are shown in the lower half of Table 6. The first and the second row in 
this table represents the results excluding and including eccentricity le into param- 
eters, respectively. Figure 18 shows that the constraints on wed become stronger by 
20% for the precessing binaries. As we stated before, imposing the prior distribution 
A/e > would increase the lower bound of wbd to the case without including /g 
into parameters. In this case, the constraint on cjbd for the precessing binaries in- 
cluding a becomes wbd > 6944, which is just 1.6 times lower than the one obtained 
ini69. 
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Fig. 18. The histograms showing the probability distributions of the lower bounds of ojbd ob- 
tained from our Monte Carlo simulations of IC* BH/NS binaries in Brans-Dicke theory^^''' . We 
take the masses of the binaries as (1.4 -|- 10^)Mq for LISA and (1.4 -|- 10)Mq for DECIGO/BBO, 
both with signal-to-noise ratio (SNR) of \/200 and we assumed that the binaries are circular. The 
thin solid and dotted lines denote the bounds obtained using LISA with and without precessions, 
while the thick dotted and solid ones represent the estimate using DECIGO/BBO with and without 
precessions, respectively. The dashed and the dotted-dashed lines represent the constraints using 
Cassini satellitel^" and PSR J1141-6545 respectively. (This figure is taken from Ref.^^^.) 

Tabic 6. The results of error estimation in Brans-Dicko theory for a (1.4-1- 1000)Mq BH/NS 
binary obtained from the Monte Carlo simulations using LISA^^''. We used two interferometers 
for the analyses and we fixed p = \/200 (p = 10 for each interferometer). The first half of the 
table shows the results for the spin-aligned case, and the second half represents the precessing 
case. The first line of each part shows the results without taking 1^ into parameters, while in 
the second line this is taken into account, a is included in all the cases. (This table is taken 
from Ref.l'^0.) 



Cases 




AlnM{%) 


A In »y 


A/3 


A In Dl 


At7s(10-^str) 


Act 


spin-aligned 
















Excluding le 


4844 


0.0396 


0.143 


0.108 


2.53 


0.406 


2.18 


Including le 


1058 


0.106 


0.241 


0.481 


1.34 


1.05 


2.45 


precessing 
















Excluding le 


6944 


0.0291 


0.107 


0.144 


0.0809 


0.341 


1.66 


Including le 


3523 


0.0432 


0.130 


0.161 


0.0851 


0.589 


1.86 



9.1.2. Proposed Constraints with DECIGO/BBO 

Next, we show the results obtained with DECIGO/BBO^''^. Here, we only consider 
circular binaries since as we explained in the previous section, the constraint on 
wbd is almost unaffected for a binary having an eccentric orbit provided that we 
impose a prior condition on eccentricity. For the instrumental noise spectrum, we 
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have used the one of DECIGO shown in Eq. (20). It is expected that BBO may give 
stronger constraints than the ones shown in this section by a factor of a few . 



Table 7. The results of error estimation with DECIGO/BBO in Brans-Dicke theory 
for various mass spin-aligned BH/NS binaries. We performed pattern-averaged esti- 
mates, using only one interferometer with p = 10. (This table is taken from Ref.^^^.) 



masses 


/in 


/fln 




AlnA^ 


Alnr; 


A/3 


Acr 




(Hz) 


(Hz) 


(106) 


(10-5) 


(%) 






(1.4-1- 1O)M0 


0.118 


100 


1.342 


0.978 


2.78 


0.190 


2.18 


(1.4-|-5O)M0 


0.0776 


85.6 


0.2662 


2.34 


2.64 


0.106 


1.09 


(1.4 + 100)Mq 


0.0651 


43.36 


0.1899 


2.34 


1.87 


0.0485 


0.563 


(1.4 + 400)Mq 


0.0460 


10.95 


0.04244 


4.96 


1.85 


0.0133 


0.250 



In Table 7, we show the pattern-averaged results of binary parameter estimation 
errors in Brans-Dicke theory with (1.4 lO)Af0, (1.4-H 5O)M0, (1.4 100)AfQ and 
(1.4 + 4OO)A/0 spin-aligned BH/NS binaries of p = 10. Just like the results obtained 
using LISA, smaller mass binaries give stronger constraints on wbd- Comparing 
these results with the ones using LISA, one sees that DECIGO/BBO has better 
ability in constraining wbd compared to LISA. This is because (i) the number of 
GW cycles are larger, (ii) the velocity at lyr before coalescence is smaller and (iii) 
the effective frequency range is larger, as shown in Table 8. 

Here, we give a rough estimate of how to derive the constraint using DE- 
CIGO/BBO. Matched filtering allows us to detect the effect due to BD correction if 
the correction term in the phase times SNR is of 0(1). From Eq. (170) and Table 8, 
we have 

^ " 84 VlO/ V5.9 X 106/ Vo.Sy V ^'lyr / V ^bd / 

(215) 

where we have used the fact that ^opn roughly corresponds to the number of 
GW cycles A^gw shown in Table 8. Therefore, if wbd is not correlated with other 
parameters, one can constrain the theory with DECIGO/BBO as Wgo'^™ > 4.3 x 10^. 
This roughly agrees with the value shown in Table 8. However, since there are 
correlations, the bound would be reduced by 100 times. This reduction is not as 
large as the one for LISA because the effective frequency range of observation is 
broader for DECIGO/BBO (as shown in Table 8). As we will discuss later, the 
expected event rate of BH/NS binaries with DECIGO/BBO would be O(104)/yr. 
This allows us to improve the bound by roughly 100 times. Therefore, in total, the 
bound would be wbd > O(IO^), which is 4 orders of magnitude stronger than the 
solar system bound 

In Table 9, we show the results of error estimation in Brans-Dicke theory for a 
(1.4 -I- 10)1^0 BH/NS binary of p = 10, for both the pattern-averaged analysis and 
of p = V200 for Monte Carlo simulations. As expected, we see that inclusion of pre- 
cession improves the constraint on wbd by a factor of two. We show the probability 
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Table 8. Comparison of the constraints on i^bd ('^bd°'^ corresponds to the bound 
assuming that cjbt> is not correlated with other parameters) and other physical 
quantities for a (1.4, 10)Mq binary with DECIGO and a (1.4, lOOO)Af0 binary with 
LISA. We performed pattern-averaged analyses and SNRs arc fixed to 10 for both 
cases. A^GW is the number of GW cycles, viy,- denotes the velocity at lyr before 
coalescence, and /in and /gn are the initial and final frequencies of observation, 
respectively. This tabic is taken from Rof.^^^. 



masses and detector 






Agw 




/in 


,ffln 




(106) 


(10«) 


(10") 




(Hz) 


(Hz) 


(1.4, 1O)M0, DECIGO 


1.34 


332 


5.9 


0.027 


0.118 


100 


(1.4,1OOO)M0, LISA 


0.00821 


21.6 


1.8 


0.083 


0.0366 


1.00 



Table 9. The results of error estimation in Brans-Dickc theory for a (1.4+10)Mq BH/NS binary 
using DECIGO/BBO. The first line shows the ones of pattern-averaged estimate. We used only 
one interferometer with p = 10. The second and the third lines show the results of Monte Carlo 
simulations. We used two interferometers for the analyses and we set p = \/200 (corresponding to 
SNR=10 for each interferometer). The second line shows the ones for the spin-aligned case, while 
the third line represents the ones including precession, cr is included in the binary parameters 
for all the cases. (This table is taken from Ref.^^^.) 



cases 




AlnM 


A In 77 


A/3 


All! Dl 


AQs 


Act 




(106) 




(%) 






(lO^^str) 




pattern-averaged 


1.342 


0.978 


2.78 


0.190 


0.100 




2.18 


no precession 


0.9774 


1.22 


3.06 


0.186 


1.24 


3.27 


2.15 


including precession 


2.317 


0.350 


0.295 


0.0551 


0.183 


2.52 


0.627 



distribution of the constraint on wbd using DECIGO/BBO in Fig. 18. The thick 
dotted histogram represents the result for the spin-aligned binaries and the thick 
soHd one shows the one for processing binaries. This shows that DECIGO/BBO 
can place 300 times stronger constraint than LISA. The reasons are the same as the 
pattern-averaged analysis. 

Unlike the case of LISA, BH/NS binaries are expected to be the definite GW 
targets for DECIGO/BBO. The event rate of BH/NS is expected to be about one 
order of magnitude smaller than the NS/NS one^^^, which is given by 10^ yr~^ 

(see also Sec. 3.3). Therefore it is possible to place even stronger constraint by 
performing a statistical analysis, where wc take the benefit of the large event rate. 
The total variance of ui is given by^^^ 

cr^2 ^ AT / ^TT[aor(z)f-h{z)—a{z)-'^dz. (216) 
Jo dz 

Here, AT = 1 yr denotes the observation time, og is the current scale factor, r{z) 
represents the comoving distance to the source, n{z) shows the BH/NS merger rate 
at redshift z and r denotes the proper look back time of the source. <j{z) represents 
the variance of Co from a BH/NS binary at redshift z. aor(z) and dr/dz are given 
in Eq. (23). n{z) is given by^^ n(z) = noR{z), where no = 10^^ Mpc^^ yr~^ is the 
estimated BH/NS merger rate today and R{z) encodes the time-evolution of this 
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rate given in Eq. (25) ^. We assume that all BH/NS binaries have the same typical 
masses of (1.4 + 1Q)Mq for simplicity. Wc first calculate the variance a{z) for each 
z using pattern-averaged estimate and obtain the total variance ctc^ using Eq. (216). 
In order to take the effects of redshift into account, we use the redshifted masses: 
"^NS ^ (1 + z)mNS and mBH -> (1 + z)m^n. 

From the pattern-averaged analysis, we find that lO"' BH/NS binaries enable 
us to place a new constraint ojbd > 2.18 x 10^ which is 94 times stronger than 
the one from a single binary of SNR=10 (or Dl = 17Gpc). This implies that if we 
perform Monte Carlo simulations, the constraint would be 

WBD > 3.77 X 10® . (217) 

This is 4 orders of magnitude stronger than the current strongest bound from the 
solar system experiments^". 

Furthermore, we have assumed that we use only two interferometers, but since 
there are eight effective interferometers in total for DECIGO/BBO, the constraint 
would improve by a factor of ■\/8/2 = 2. Moreover, we note again that if we use 
BBO which is roughly three times more sensitive than the noise curve that we used 
in this section, the constraint would further be three times stronger. In total, BBO 
should be able place nearly 5 orders of magnitude stronger constraint than the solar 
system bound! This is still 2 orders of magnitude stronger than the expected bound 
from the future solar system mission ASTROD^^^. 



9.1.3. Prospects for Other Space- Borne GW Interferometers 

In this subsection, we explain the proposed bounds on wed using eLISA and 
ASTROD-GW. The non-sky-averaged instrumental noise spectral density for eLISA 
is given by®^ 



-iinst,oLISA 



(/) = 



4^acc(/)/(27r/)4 + SUf) + ^om„(/) 



LI 



1 + 



f 



0.41c/(2Le) 



(218) 



^R{z) is taken from the redshift evolution of NS/NS merger rate obtained in Ref.^^^. Although 
BH/NS merger rate has not been obtained in this referenee, the formation rate of BH/NS follows 
the same evolution as that of NS/NS (see Fig. 6 of Ref.-^-'^). Therefore wc adopt this evolution for 
BH/NS merger rate as well. 

'ET is expected to detect BH/NS signals out to z ^ 1 for the SNR threshold value of 10 and 
the expected event rate is more than O(10'')/yr. Following the analysis of DECIGO/BBO, we 
found that ET can place ljbd > 7.25 x 10^ using a single interferometer for the pattern-averaged 
analysis. (See Ref. ^-^^ for the proposed constraint on ojbd with ET without including spins into 
parameters.) If the event rate is about 10 times larger (which corresponds to the averaged value 
found in e.g. Ref.^^) and if we assume that we use 2 interferometers, the above constraint further 
improves as lobd > 3.24 X 10^. This is roughly 1 order of magnitude stronger than the current 
solar system bound. For processing binaries, this bound should improve by a factor of a few. 



12, 2013 2:23 



WSPC/INSTRUCTION FILE papcr-arXiv 



60 Kent Yagi 



Table 10. The proposed bounds on ojbd with 
LISA, eLISA and ASTROD-GW. Wc assume 
spin-aUgned BH/NS binaries of p = 10 with a single 
interferometer and perform pattern-averaged Fisher 
analysis. We set NS masses to 1.4Mq and include 
both (9 and a into parameters, while we consider 
binaries with circular orbits and do not include le 
into parameters. The LISA bounds are same as the 
ones shown in Table 5. 



BH mass 


LISA 


cLISA 


ASTROD-GW 


400 Mq 


24886 


25644 


24809 


1000 Mq 


8210 


8759 


8166 


5000 Mq 


1933 


2520 


1861 


10000 Mq 


862 


1148 


820 



where = 1 Mkm the acceleration noise, shot noise and all other measurement 
noises (for the position noises) are given by 

5acc(/) - 2.13 X 10-29 (^1 + Hz"! , (219) 

Ssnif) ^ 5.25 X 10-23 Rz"^ , (220) 
^omn(/) = 6.28 X 10-23 ni^ R^-' , (221) 

respectively. For ASTROD-GW, the non-sky-averaged instrumental noise spectral 
density is given by^^'' 

^inst,ASTROD ^ [1 + jf / fafmSp + 4SJ [2^ f)^ ^ (222) 

where fa = c/(27rLa) with La = 52 Mkm and the acceleration and the position 
noises are given by 

Sa = 1.0816 X 10-1* m^ Hz-i , (223) 
Sp = 9.0 X 10-^" m2 Hz-i , (224) 

respectively. The total noise sensitivities including the WD/ WD confusion noises 
can be obtained by substituting Eqs. (218) or (222) into Eq. (18), and they are 
shown in Fig. 2. 

The bounds using these detectors, together with the ones with LISA, are shown 
in Table 10. One can sec that the proposed constraints arc almost the same among 
these detectors. (Ref.-'^^s discusses that the bound from eLISA is worse by a factor 
of a few compared to the LISA one. This difference comes from the difference in 
the observational frequency range that we use for eLISA.) Since we are setting SNR 
to be p = 10, it is the shape of the sensitivity curves that matters and not the 
overall amplitude of the curves. For the mass range of BHs shown in Table 10, the 
frequency at lyr before coalescence is higher than lO-^Hz. For this frequency range, 
the shape of the sensitivity curves of the 3 detectors arc similar (a/S^ (x f^), and 
hence the bounds turn out to be almost identical. The bounds would be similar even 
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if we include the effect of precessions, as can be understood by comparing Tables 5 
and 6. 

9.2. Massive Gravity 

9.2.1. Proposed Constraints with ground-based detectors, LISA, Pulsar 
Timings and CMB 

In massive gravity theories, the phase velocity of GW is given as"'^*^ 



We see that it depends on its frequency, which modifies the time of arrival from 
general relativity. This modifies the phasing of the gravitational waveforms. Jones^^"'^ 
estimated how accurately LISA will be able to distinguish the difference in the 
arrival times of different harmonic signals from a binary with an eccentric orbit. 
Larson and Hiscock^^^ proposed to place bounds on Xg by using both GW and EM 
signals from a WD binary. This work was improved by Cutler et al.^^^ and Cooray 
and Seto^^'^ where the latter considered ~400 close WD binaries and obtained the 
bound Ag > 2.1 X lO^^cm which is two orders of magnitude stronger than the solar 
system experiment. Kocsis et al^^^ proposed to use the correlation between GW and 
EM signals of SMBH binary coalescence. Adopting the timing uncertainty as the 
inverse of the GW frequency at ISCO, they estimated the possible future bound as 
Ag > 2.8 X lO'^^cm, which is three orders of magnitude stronger than solar system 
bound. However, uncertainty remains on the systematic delay in the emission of 
EM bursts. They claim that if the variability of EM signals before coalescence is 
identified and can be related to the orbital period, similar analyses as the ones 
using WD binaries mentioned above can be applied so that the systematic errors 
arc reduced and one can obtain stronger constraint. These results are summarized 
in Table 11. 




(225) 



Table 11. The proposed bounds on Ag 
using GWs. 



Current Bound 



Ag(cm) 



EM counterparts 
Pulsar Timing 



Eccentric Binary 
WD 




Will^®'' included Ag into the binary parameters and performed the matched fil- 
tering analysis, estimating how accurately one can determine Ag using adv. LIGO 
or LISA. Will and Yunes^^® followed similar analysis using the improved noise curve 
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for LISA. As in the BD case, they did not include the spins and also they only per- 
formed the pattern-averaged analysis. Berti et al}^^ estimated the constraint on Xg 
by performing the Monte Carlo simulations, taking the spin-orbit coupling into ac- 
count. They find that for a (10*^ + 10*^)Mq BH/BH binary at 3Gpc, 1 yr observation 
with LISA can place the bound Xg > 1.33 x lO^^cm on average. Arun and Will^^® 
included the effect of higher harmonics and found that for higher mass sources, the 
constraints become stronger compared to the analysis using the restricted wave- 
forms. The above results are summarized in the upper halves of Tables 12 and 13. 

Table 12. Summary of the works on probing MG theories using GWs from compact binaries 
with the matched filtering analysis. The second column denotes up to which PN order they take. 
(Berti et al. (2011) is left blank because they do not state this information in their paper.) 
The third, fourth, fifth and sixth columns represent whether they take the spin-orbit coupling /3, 
higher harmonics (HH), precession or eccentricity into account. The seventh column shows whether 
they consider all of the inspiral, merger and ringdown (IMR) phases or just the inspiral phase. 
The eighth column describes whether they assume multiple-source detections. The ninth column 
represent whether they perform model-independent (MI) calculations. Finally, the last column 
shows what kind of analyses they have performed (either the pattern-averaged (PA), Monte Carlo 
(MC) or Bayesian (B).) 



Reference 


PN 


/3 


HH 


prec. 


ecc. 


IMR 


multi. 


MI 


analy. 


Will (1998)^«''' 


1.5 


X 


X 


X 


X 


X 


X 


X 


PA 


Will & Yuncs (2004)188 


1.5 


X 


X 


X 


X 


X 


X 


X 


PA 


Berti et al. (2005) ^''^ 


2 


o 


X 


X 


X 


X 


X 


X 


MC 


Arun & Will (2009) 


3.5 


o 


o 


X 


X 


X 


X 


X 


PA 


Stavridis & Will (2009)19° 


2 


o 


X 


o 


X 


X 


X 


X 


MC 


Yagi & Tanaka (2010)"° 


2 


o 


X 


o 


o 


X 


X 


X 


MC 


Yagi & Tanaka (2011)"^ 


2 


o 


X 


o 


o 


X 


X 


X 


MC 


Keppel & Ajith (2010)i«9 


3.5 


X 


X 


X 


X 


o 


X 


X 


PA 


Del Pozzo et al. (2011)1^1 


2 


X 


X 


X 


X 


X 


o 


X 


B 


Cornish et al. (2011)1*5 


3.5 


X 


X 


X 


X 


X 


X 


o 


B 


Berti et al. (2011)1^2 




X 


X 


X 


X 


X 


o 


X 


MC 


Huwyler et al. (2011)14'^ 


2 


o 


o 


o 


X 


X 


X 


o 


MC 



Following the BD case, we improved the analysis of Ref.^^^ by (i) including 
the spin-spin coupling, (ii) considering binaries with slight eccentricity and (iii) 
looking at precessing binaries. In Table 14 we show the pattern-averaged results for 
various spin-aligned BH/BH binaries. One sees that including both a and le into 
parameters increases the measurement error of Xg only slightly, implying that /3g is 
not so strongly correlated with other parameters. In the massive gravity case, Xg is 
more degenerate with a than with ■ Unlike in the Brans-Dicke theory, both terms 
containing Xg and a in the phase ^'(/) have frequency dependences with positive 
power-law indices. Therefore, their correlation is stronger than the one between Xg 
and le- Notice that binaries with larger masses place stronger constraint on Xg. This 
is because the correction to Uph is larger when the frequency is lower. 

Next, we show in Table 15 the results of Monte Carlo simulations. The upper 
half shows the ones for the spin-aligned binaries. We fix the binary mass to (10^ + 
10^)Mq, and the distance to Dl = 3Gpc. The corresponding histogram (only for 
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Table 13. Summary of the constraints on Ag using GWs. The units in the table are lO^^cm for adv. LIGO 
and ET, while lO^^cm for LISA, eLISA and DECIGO/BBO. The numbers in brackets denote the total 
masses of the binaries in the unit of A/q . 



Reference 


adv. LIGO 


ET 


LISA 


cLISA DECIGO/BBO 


Will (1998)^**'' 


0.6 (20) 




7 (2 X 10 




Will & Yunes (2004)168 






5 (lO'^) 




Berti et al. (2005) I'^S 






1 (2 X lO'') 




Arun & Will (2009)1*8 


0.7 (60) 


10 (400) 


5 (2 X 10*^) 




Stavridis & Will (2009)1^0 






7 (2 X lO'') 




Yagi & Tanaka (2010) i^° 






3 (1.1 X lO'^) 




Yagi & Tanaka (2011)1^^ 








0.3 (1.1 X lO'') 


Kcppel & Ajith (2010)1**^ 


8 (360) 


70 (3000) 


60 (4.8 X 10^) 




Del Pozzo et al. (2011)i''l 


0.5-2.5 








Cornish et al. (2011)"5 


0.9 (18-24) 




4 (4 - 5 X 10^) 




Berti et al. (2011)1^2 






6.5-7.5 


3-5 


Huwyler et al. (2011)1*'' 






7 (1.3 X 10'') 





Table 14. The results of error estimation in massive gravity theories using LISAl™. These 
are calculated with pattern-averaged analysis for spin-aligned SMBH binaries with masses 
(10'^ -I- 1O'')M0, (10'' -I- 10^)Mq, (id" -I- 106)Mq and (10*5 -I- 10^)Mq at 3Gpc. We used on 
only one interferometer. As in the Brans-Dicke case, the first line of each binary represents the 
estimation without taking cr nor /g into parameters. These results do not exactly match with 
the ones shown inl^^ because we also included prior information. The meaning of second and 
third lines arc the same as in Table 5. (This table is taken from Rcf.i'"''.) 



(J 


le 


As 

(10^" cm) 


AlnM 

(%) 


A buy 


A/3 


Ate 

(s) 


A<Pc 


Acr 


Ale 

(10-10) 


(10' 4- 

X 


10'')Mq 

X 


22.77 


0.0669 


0.467 


2.93 


75.7 


1.06 






o 


X 


11.33 


0.0687 


0.960 


7.10 


77.8 


1.09 


1.73 




o 


O 


11.29 


0.246 


1.10 


7.56 


133 


2.23 


1.89 


13.0 


(10 ' 4- 

X 


1O'')M0 

X 


9.629 


0.0493 


0.253 


1.49 


82.3 


2.15 






o 


X 


4.061 


0.0495 


0.839 


6.26 


82.8 


2.16 


1.89 




o 


o 


4.052 


0.202 


0.927 


6.54 


159 


5.01 


1.94 


14.3 


(10'' 4- 

X 


10'')Mo 

X 


12.41 


0.00869 


0.122 


0.787 


3.01 


0.316 






o 


X 


3.582 


0.00871 


0.926 


6.97 


3.01 


0.316 


1.68 




o 


o 


3.582 


0.0288 


0.933 


7.00 


4.75 


0.589 


1.69 


1.14 


(10« 4- 

X 


1O'')M0 

X 


6.019 


0.00586 


0.0551 


0.337 


2.45 


0.521 






o 


X 


1.286 


0.00586 


0.823 


6.20 


2.45 


0.521 


1.88 




o 


o 


1.285 


0.0200 


0.826 


6.21 


3.96 


0.989 


1.88 


1.50 



circular binaries) for the probability distribution of Ag with LISA is shown as the 
(light blue) thin dotted one in Fig. 19. When including higher FN terms, the Fourier 
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Table 15. The results of the Monte Carlo simulations in massive gravity theories for (10^ + 10^)Mq 
BH/BH binaries at 3Gpc without pattern averaging using LISA^^". We used two interferometers for the 
analyses. As in Brans-Dicke case, we distribute 10^ binaries, calculate the error of each parameter for 
each binary and take the average. The first half of the table shows the results for the spin-aligned cases, 
and the second half represents the ones for the processing binaries. The first line of each part shows the 
results without taking le into parameters, while in the second line this is taken into account, a is included 
in the parameters for all the cases. (This table is taken from Ref.^^''.) 



Cases 


As 


SNR 


AlnM 


Alnri 


A/3 


AlnDL 


AQs 


Act 




(lO^lcm) 




(%) 






(%) 


(lO^-^str) 




Spin-Aligned 


















Excluding le 


0.40598 


1540 


0.0507 


0.841 


6.27 


2.30 


0.957 


1.89 


Including le 


0.40507 


1540 


0.191 


0.927 


6.54 


2.33 


0.972 


1.94 


Precessing 


















Excluding le 


4.8540 


1596 


0.00838 


0.00675 


0.0117 


0.189 


0.366 


0.0508 


Including le 


3.0570 


1586 


0.0269 


0.00708 


0.0120 


0.192 


0.364 


0.0825 




Fig. 19. The histograms showing the probability distributions of the lower bounds of Xg obtained 
from our Monte Carlo simulations of 10* BH/BH binaries in massive gravity theories^'^"' ^'^^ . We 
take the masses of the binaries as (lO'^ + 10'^) Mq for LISA and {10^ + 10^)Mq for DECIGO/BBO, 
both placed at the luminosity distance of 3Gpc. The meaning of each histogram is the same as in 
Fig. 18. (This figure is taken from Ref.-'''^.) 



transform of the GW waveform for a circular binary can be expressed as 



*(/) = 27r/io -cl)o-- + *opn(./) 



, 128^ n/. /3715 55 



-4(47r - P)x^/^ + 



/ 15293365 27145 



3085 



77^ - lOcr ) .T^ 



(226) 



V 508032 504 ' ' 72 
with X = {jTAItfY^'^ ■ Ahhough there are 5 parameters {Ai,r], I3,a, and /3g) inside 
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the square brackets, there are only 4 PN terms in total (the leading quadrupole 
term, IPN, 1.5PN, and 2PN), meaning that these parameters degenerate for the 
spin-aligned binaries. Since the level of uncertainty of a is determined completely 
from the prior distribution, the uncertainty in Pg is also determined by the one in 
(T, and hence there is a sharp peak in the histogram in Fig. 19. 

On the other hand, the lower bounds on Ag for precessing binaries are shown 
in the second half of Table 15. We fix the masses of the binary component to 
(10^ + 1O^)M0 and assume the BH spins to be 0.5 and 0. Since one of the BH spins 
is 0, the situation reduces to the simple precession. The corresponding histogram is 
shown as the (blue) thin solid one in Figs. 19. One sees that the lower bounds on 
Ag increase by one order of magnitude for the precessing binaries. The lower tail in 
the histogram corresponds to the binaries with n 1, i.e. the binaries are almost 
spin-aligned. 

Notice the impact that the precession has on massive gravity theories compared 
to Brans-Dicke theory. As explained above, up to 2PN order, there are not enough 
number of PN terms in GW phase to determine all of binary parameters. However, 
the spin precession brings additional information that helps in determining the 
spins, which solves the degeneracies between spins and other parameters. Therefore, 
inclusion of spin precession is crucial when constraining massive gravity theories. 
On the other hand, there are enough number of PN terms contained in Brans-Dicke 
theory, so that the effect of the precession is relatively weak. 

Just before wc completed our paper^™, similar paper by Stavridis and Will^^'' 
appeared. Although we assumed the simple precession which can be applied only to 
binaries with equal-mass components or one of the spins are zero, they numerically 
solved precession equations so that their estimate is not limited to such binaries. 
For precessing binaries of (10^ -I- 1O'')M0 at 3Gpc, LISA can place the bound Ag > 
5 X lO^^cm which is slightly larger than our estimate. They found an interesting 
result that this constraint is almost the same as the one without spins, suggesting 
that precession makes the spin and other parameters degenerate. Very recently, 
calculations including both precessions and higher harmonics has been performed 
by Huwyler et al}'^'' . Their constraint on Ag is comparable to the one by Stavridis 
and Will. It seems that the precessions already solves the degeneracies between Ag 
and other binary parameters so that the effect of higher harmonics is not so strong. 

Keppel and Ajith^®^ calculated the bound including merger and ringdown in- 
formation for the first time. They used the phenomenological inspiral-merger- 
ringdown waveform^°^'^°^ (see Sec. 5.3). They performed the pattern-averaged anal- 
ysis and obtained rather stringent bound Ag > 5.9 x lO^^cm for the total mass of 
Mt = 4.8 X 10'' M© equal-mass binary at 3Gpc with LISA, and Ag > 7.8 x lO^^cm 
for the total mass of Mt = 36OM0 equal-mass binary at IGpc with adv. LIGO. The 
latter result shows that the second-generation ground-based interferometers can put 
the bound a few times stronger than the current solar system one. Although their 
analysis does not include spins, calculations by Stavridis and Will implies that the 
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results should not change mueh when one includes spins and precession effects. 

Instead of using Fisher analyses, Del Pozzo et al}^^ performed Bayesian analyses 
and found that the second-generation ground-based interferometers can put the 
bound Xg > 2.6 x lO^^cm using 50 inspiral signals. Similar Bayesian analysis has 
been performed by Cornish et al.^'^^ for model-independent test of modified gravities 
applying PFE formalism. Their constraint is consistent with the one found by Del 
Pozzo et al}^^ . 

Most of the works mentioned above assumes a binary of masses Mt ~ lO^M© 
at distances ^ 3Gpc for the constraint with LISA. However, BH merger tree sim- 
ulations show that the typical binaries that are observable with LISA have smaller 
masses with larger distances^®^"^^"'^ . Recently, Berti et al}^^ performed 1000 real- 
izations of merger simulations for two BH formation scenarios, and found that for 
each realization, there are about 40 BH binary mergers that can be observed with 
LISA. By combining signals from these detectable binaries, they found the con- 
straint Xg > 6.5-7.5 X lO^^cm with LISA. They also estimated the bound using 
newly-proposed LISA (cLISA) and found that the bounds are reduced roughly by a 
factor of two. These results are summarized in the lower halves of Tables 12 and 13. 

Apart from GW observations using interferometers, there are works that pro- 
pose bounds on Xg using pulsar timings. Baskaran et al.^^^ considered the so-called 
surfing effect"^^"^ which is the resonance effect arising when the speeds of light and 
GWs are unequal. From current pulsar timing observation of PRS B1937-I-21 with 
the observation period Tobs = 4yr and the timing residual i?rms = 0.1//sec, they ob- 
tained the bound Xg > 1.5 x lO^^cm, assuming that the energy density of the SMBH 
stochastic GW background is^^^ f^cwCT^bs) = 4.2 x 10"^" and the spectral tih^^^ 
riT = 2/3. Due to this assumption, this bound is less robust compared to the solar 
system experiment. By assuming the future 10-yr pulsar timing arrays of 300 pulsars 
with 100ns timing accuracy, Lee et al^^^ estimated the effect of massive graviton 
on the correlation between pulsar signals and obtained Xg > 4.1 x lO^^cm. Also, the 
graviton mass can be investigated by CMB observations. For certain ranges of Xg, 
the graviton mass has a remarkable effect on the large scale GMB power spectrum 
(for both temperature fluctuations and polarizations) The results mentioned 

in this paragraph are also included in Table 11. 

9.2.2. Proposed Constraints with DECIGO/BBO 

In this section, we show the results using DECIGO/BBO. Following Sec. 9.1.2, we 
have adopted DECIGO instrumental noise rather than the BBO one. In Table 16, 
we show the pattern-averaged results for various spin-aligned BH/BH binaries. As 
discussed in the LISA section, since the constraints would be stronger for larger 
mass binaries, LISA performs better than DECIGO in constraining massive gravity 
theories. (The reason a (10® -f lO^)A/0 binary gives stronger constraint than a 
(10® -I- lO®)Af0 binary is because the GW signal shifts out of the observation band 
as we increase the mass.) 
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Table 16. The results of error estimation in massive gravity theories for spin-aligned BH/BH bina- 
ries at 3 Gpc with various masses using DECIGO/BBO. We performed pattern- averaged analyses 
using only one interferometer. (This table is taken from Ref.^^^.) 



masses 


/in 


/fin 


SNR 


As 


AlnX 


Aln»y 


A/3 


Act 






(mHz) 


(mHz) 




(102"cm) 


(%) 








(10'' 


+ 1O'')M0 


1.0 


2.20 


1338 


1.014 


14.0 


2.46 


9.40 


2.50 


(106 


+ 1O5)M0 


1.0 


4.00 


2044 


1.270 


1.19 


1.69 


9.10 


2.44 


(105 


+ 1O5)M0 


1.0 


22.0 


4909 


1.133 


0.0286 


0.930 


7.00 


1.69 


(105 


+ 10'i)Mo 


1.0 


40.0 


3021 


0.4066 


4.51x10"^ 


0.823 


6.20 


1.88 




+ 10-4)1/0 


1.0 


220.0 


29569 


0.3852 


3.54 xlO"* 


0.924 


6.96 


1.68 



Table 17. The results of error estimation in massive gravity theories for (10® -|- lO^)Af0 
BH/BH binaries at 3Gpc using DECIGO/BBO. The meaning of each line is the same as in 
Table 9. (This tabic is taken from Ref.^'^S ) 



cases 


SNR 


As 


AlnA^ 


Alnr? 


A/3 


Ao- 


A^s 






(1020cm) 


(%) 








(str) 


pattern-averaged 


2044 


1.270 


1.19 


1.69 


9.10 


2.44 




no precession 


2601 


1.266 


1.16 


1.64 


8.92 


2.40 


1.16 


including precession 


2666 


3.349 


0.314 


0.0388 


0.0612 


0.529 


0.0248 



In Table 17, we show the results of binary parameter estimation in massive 
gravity theories for both pattern-averaged analysis and Monte Carlo simulations 
with and without including spin precession. We assume (10^ + 10^)A'/o BH/BH 
binaries at Dl = 3Gpc. The meaning of each row is the same as in Table 9. One sees 
that the constraint on Xg becomes twice stronger when we include precession. The 
corresponding histograms are shown in Fig. 19. Notice that the effect of precession 
is larger for LISA. We expect that this is due to the wider effective frequency range. 
The proposed lower bound on Xg with DECIGO/BBO is Xg > 3.35 x lO^Ocm. 
Although this is one order of magnitude weaker than the LISA one, it is still three 
orders of magnitude stronger than the current solar system bound^®^. 

9.2.3. Proposed Constraints with DPF 

In this section, we show the results obtained using DPF. First, we give a rough 
estimate on how strong we can constrain Xg with DPF. It is not possible to detect 
the effect of massive gravity if the correction term in Eq. (186) is smaller than p~^. 
This leads to the constraint 

^ - \5kpcJ Vo.iHzy V 30 y ^ ^ 

Next, we estimate the constraint numerically using Fisher analysis for a spin- 
aligned binary using inspiral-mcrgcr-ringdown hybrid waveform. In Fig. 20, we show 
the sky-averaged results for the lower bound on Xg against various BH masses, 
assuming spin-aligned, equal-mass binaries with circular orbits, with x ~ 0-2 at 
Dl = 5kpc. For a (10^ + lO'')A/0 binary, we find the constraint of Ag > 6 x lO^^cm 
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DPF, withx 

adv. DPF, with % 

DPF, without X 

adv. DPF, without x 

sol. sys. exp. 
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10"^ 10^ 
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Fig. 20. The lower bounds on Aj using DPF with (red thick soUd) and without (magenta thick 
dotted) taking x into parameters, and using adv. DPF with (green thin solid) and without (blue 
thin dotted) taking x into parameters. We assumed spin-aligned, equal-mass binaries with x = 0-2 
at Dl = 5kpc. We also assumed that the orbits are circular. The (black) thin dotted horizontal 
line at Ag = 2.8 X lO^'^cm represents the bound obtained from solar system experiments*^ while 
the (black) thick dotted-dashed horizontal line at Ag = 1.6 X lO^^cm shows the one obtained from 
binary pulsar test^^'^. (This figure is taken from Ref. s^*.) 

using DPF, which is weaker than our rough estimate in Eq. (227). This is due to 
the degeneracies between f3g and other binary parameters. In Fig. 20, we also show 
the constraints without including x into variable parameters (i.e. assuming that x 
is know a priori). This roughly gives how well one can constrain the theory when 
the binary is precessing^^^. This corresponds to the case where the degeneracies 
between Ag and x are completely disentangled. The dotted-dashed horizontal line 
at Xg ~ 2.8 X lO^^cm represents the (static) lower bound obtained from the solar 
system experiments®^. Although DPF constraint is slightly weaker, they are still 
meaningful since DPF measures the deviation in the propagation speed of GWs that 
appear in gravitational waveform phase, while the solar system experiment measures 
the deviation in the effective gravitational constant (or in the Kepler's third law). 
We find that DPF can place about 2 orders of magnitude stronger (dynamical) 
constraint than the one from the binary pulsar^'''^ in the weak-field regime, shown 
as the dotted-dashed horizontal line at Ag = 1.6 x lO^^cm. Furthermore, Finn and 
Sutton^''^ assumed Fierz-Pauli-type theory while the constraint obtained here is 
independent of the specific type of massive gravity theories. Since the effect of the 
finite mass of graviton is larger for larger mass binaries, it would be important to 
reduce the acceleration noises, just like in the case of joint search with DPF and 
the network of ground-based detectors. 

In Fig. 21, we show the probability distribution of the lower bound on Ag with 
Monte Carlo simulations. We assume 10" Af© equal -mass binaries of x — 0-2 with 
circular orbits in co Centauri. One sees that the constraint is now slightly weaker 
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Fig. 21. Histograms showing the lower bounds on Ag obtained from Monte Carlo simulations 
with 10* M0 equal-mass IMBH binaries with circular orbits of x = 0.2 in u) Centauri with random 
orientation for the orbital angular momenta. We show here the results using DPF. The (black) 
thick and (red) thin histogram each represents the one with and without taking x as a variable 
parameter, respectively. (This figure is taken from Ref.^^*.) 

than the one with the sky-averaged analysis. However, this is still much stronger 
than the one from the binary pulsar observations. These probability distributions 
shift to larger Ag by roughly y/l = 2.6 when we use adv. DPF. 

9.2.4. Prospects for Other Space- Borne GW Interferometers 

Tabic 18. The proposed bounds on Ag (in units of 
10^1 cm) with LISA, eLISA and ASTROD-GW. We as- 
sume spin-aligned BH/BH binaries at 3Gpc with a sin- 
gle interferometer and perform pattern-averaged Fisher 
analysis. We exclude both /3 and a into parameters. 
(This should mimic the precessing case.) We consider 
binaries with circular orbits and do not include 7e into 
parameters. 



BH 


masses 


LISA 


cLISA 


ASTROD-GW 


(10 ' -t 


- W)Mq 


6.941 


1.914 


12.79 


(loM 


- 106)Mq 


3.946 


1.300 


5.851 


(10« 4 


- 1O6)M0 


4.631 


2.010 


5.411 


(10*5 4 


~ 1O5)M0 


2.743 


1.330 


2.492 



In this subsection, we explain the proposed bounds on Ag using eLISA and 
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ASTROD-GW. The bounds using these detectors are shown in Table 18. As for the 
BD case, we show the resuhs with LISA for comparison. Here, we do not include 
spins into parameters. The bounds for the precessing case should be similar, as can 
be understood by comparing Tables 15 and 18. eLISA can place the bounds that 
are weaker than the LISA ones by a factor of a few. This is because eLISA is less 
sensitive in the lower frequency part. Among the GW interferometers considered 
in this paper, ASTROD-GW performs the best in constraining massive gravity 
theories. This is because it is the most sensitive in the lower frequency range. It can 
place even stronger constraint if it detects GW signals from a binary with larger 
BH masses. 

9.3. Other Theories 

Up to here, we focused on BD and MG theories, but GW observations will allow 
us to test other theories. One example is quadratic gravity^^^ in which general 
quadratic curvature terms coupled to a scalar field is added to Einstein-Hilbert 
term at the level of action, together with the kinetic and the potential terms of the 
scalar field. This theory includes Einstein-Dilaton-Gauss-Bonnet (EDGB)^^'' theory 
in the even parity sector and dynamical Ghern-Simons (CS)^^^''^^^ gravity in the 
odd parity sector as specific examples, which are motivated from e.g. superstring 
theory-'^'^'^' The gravitational and scalar radiation energy flux in quadratic gravity 
are calculated in Ref.^^*^ within the post-Newtonian approach. In the even parity 
sector, the scalar dipole radiation exists, just like in BD theory, and this gives -IPN 
dissipative correction to the gravitational waveform phase relative to GR. This 
dominates the OPN conservative correction coming from the scalar interaction. In 
Ref.'^'^^, the current strongest bound on EDGB theory is obtained from the orbital 
decay rate of the low-mass X-ray binary A0620-00'^''^, and possible future constraints 
using GW interferometers are also estimated. Ground-based detectors are not so 
useful in constraining the theory, but eLISA should be able to place comparable 
bound and DEGIGO/BBO can place 100 times stronger constraint compared to 
the current bound. 

As for the odd parity sector of quadratic gravity (dynamical GS gravity), both 
scalar radiation and corrections to the gravitational radiation give 2PN dissipative 
correction to the gravitational waveform phase'^™. Also in this theory, there is a 
quadrupole moment deformation in the BH solution valid to quadratic order in 
spin,^*^^ which gives 2PN conservative correction. Another conservative correction 
at the same PN order also arises from the scalar interaction. Gombining all of 
these 2PN CS effects, we construct gravitational waveform from BH binaries and 
performed Fisher analysis^*^^ to obtain the constraints on the characteristic length 
scale of the theory (The theory reduces to GR in the limit ^ 0.) We present 
proposed constraints on the theory with GW observations in Fig. 22. We found that 
ground-based detectors and DECIGO/BBO have ability to constrain the theory 
by 6-7 orders of magnitude stronger than the solar system bound'^°^ of ^^/^ < 
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Fig. 22. Projected Icr upper bounds on the characteristic length scale of dynamical CS gravity 
^'^Z* in km with the second-generation ground-based detectors (top left), ET (top right), DE- 
CIGO/BBO (bottom left) and LISA (bottom right) for spinning BH binaries with mi/m2 = 2. 
The fixed values for the dimensionless spin parameter of the primary xi and the luminosity dis- 
tance Dij are shown at the top right of each panel. The constraints using the second-generation 
ground-based detectors have been obtained by assuming that the spins are known a priori. This 
roughly models projected constraints with precessing BH binary observations and may be cor- 
rect within an error of 30%. For other detectors, binaries are assumed to be spin-aligned (or 
anti-aligned). The colored contours show the regions of parameter space where the constraints 
on ^1/'* also satisfy the small coupling approximation (the approximation that we used to derive 
the gravitaional waveform) at the fiducial luminosity distances chosen. This figure is taken from 
Ref.303. 



C'(10^)km. On the other hand, LISA is not as good as these detectors, but stiU, it 
should be able to place 2-3 orders of magnitude stronger constraint than the solar 
system experiments, which is roughly consistent with the resuhs in Refs.'^"^''^^^. 
(The difference comes from the fact that Ref.'^^'^ considers 2PN dissipativc and 
conservative effects at quadratic order in spin while Rcfs.'^^''^'^'^® consider 4.5PN 
conservative effect at linear order in spin.) If the BHs arc not spinning, then the 
corrections to the scalar and gravitational radiation flux would be of 7PN and 
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Fig. 23. The scalar radiation energy flux E^^^ relative to the GR gravitational one Egk against 
v/c for a non-spinning BH binary in dynamical CS gravity. The crosses are the ones obtained 
numerically in Ref.^'^* using BH perturbation while the solid curve is the analytic result obtained 
in Ref.^"" using PN formalism. Observe that at low velocities, in the regime where the PN ap- 
proximation is valid, the two curves agree. (This figure is taken from Ref.^"''.) 



6PN orders respectively. This is first calculated numerically by Pani et al. using 
BH perturbation. Later, the scalar radiation correction is confirmed analytically 
under PN analysis (as shown in Fig. 23) and the origin of 6PN correction for the 
gravitational radiation is also discussed^°°. Quasi- normal modes in this theory for 
non-spinning BHs have been calculated by Molina et al.^'^^ and the corrections might 
be detected by future GW observations. 

For other theories, gravitational waves from compact binaries allow us to probe 
the size of the extra dimension'^^^'^^^ and the time variation of the gravitational 
constant G GWs from NS oscillations can be used to constrain not only scalar- 
tensor theories as mentioned above, but also other theories such as tensor-vector- 
scalar theory^-'^'^"^-'^^. Direct detections of additional polarizations^-'^^''^^^ would be 
smoking-guns for the deviation from GR. 

10. Conclusions 

Observations complementary to the ground-based GW interferometers can be per- 
formed with the space-borne ones. As for the prototype of DECIGO, DPF is aimed 
to be launched in 2016-2017. The main goals of this satellite are to test the im- 
portant technologies that are crucial to the space GW mission, and to carry out 
observations of GWs and Earth gravity. Although the event rate is not so promis- 
ing, DPF is sensitive enough to detect GWs from galactic IMBH binaries. In this 
review article, we first showed what information can be achieved if these GWs are 
detected. Since the sensitivity frequency range of DPF is lower than the ground- 
based detectors, DPF can detect inspiral signals earlier than the latter. Therefore, 
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DPF may be able to give an alert to the ground-based detectors so that they can 
be prepared for observing merger and ringdown phases. Also, DPF can determine 
binary parameters such as masses and spins to several % accuracies. On the other 
hand, the measurement accuracy of the distance and the angular resolution arc not 
so high. 

One very important aspect of GW observations is that they can test alternative 
theories of gravity in the strong-field regime. In this article, we considered extending 
2 of the most important characteristics in GR: (i) there are only tensor degrees of 
freedom and (ii) the graviton is massless (or GWs propagate at the speed of light). 
We consider BD and MG theories as representatives of the modifications of (i) and 
(ii), respectively. The former is the simplest example of the scalar-tensor theory. 
There exists a scalar dipole radiation which changes the evolution of the binary 
system. This gives a correction to the GW phase at "-1PN" order relative to GR. It 
has a negative PN order correction because scalar radiation is dipolar while GWs are 
quadrupolar. This correction is suppressed by the inverse of the BD parameter that 
describes the strength of the coupling between the scalar and the matter fields. We 
carried out the Fisher analysis and found that DECIGO/BBO can place 4-5 orders 
of magnitude stronger constraint compared to the solar system experiment. This 
constraint is still 1-2 orders of magnitude more stringent than the one from future 
solar system mission ASTROD^''^ DECIGO/BBO has more advantages than LISA 
because (i) the number of GW cycles are larger which allows us to perform more 
accurate test, and (ii) the velocity of the binary constituents is smaller which gives 
larger dipolar correction. We also found that eLISA and ASTROD-GW would place 
similar bounds as LISA if we fix the SNRs to be the same. We note that scalar-tensor 
theories might also be constrained from GWs associated with NS oscillations'^^^''^^^. 

As for MG theories, the propagation speed of GWs now becomes smaller than 
the speed of light, which gives IPN correction to GW phase. We carried out similar 
analysis mentioned above and found that in this case, LISA performs better than 
DECIGO/BBO. This is because the correction to the propagation speed becomes 
larger for lower frequency GWs. LISA can place 4 orders of magnitude stronger con- 
straint than the solar system bound. DPF may be able to place comparable bound 
to the solar system one and 2 orders of magnitude stronger constraint than the 
binary pulsar observations. Space-borne GW interferometers would supply valuable 
information on astrophysics and gravity. eLISA would place slightly weaker con- 
straint than LISA, while ASTROD-GW can place the strongest bounds among all 
the GW interferometers considered in this paper. 

One might worry that the magnification effect of the gravitational lensing on 
GWs'^^^''^^-'^ might spoil the results shown in this article. However, this effect only 
affects the amplitude of GWs, and hence, as discussed e.g. in Ref.^'*, this increases 
the measurement errors on the luminosity distance and the angular resolution of the 
binary. However, since the amplitude and the phase are almost uncorrelated, this 
effect would not affect the measurement accuracies of the parameters that appear 
in the phase. 
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Appendix A. Lengthy Formulas 
Appendix A.l. Spin-Aligned Binaries 

When expressing the observed waveform for the spin-aligned case, we need to express 
9s{t), 0s (i) and ipsit) in terms of the angles ^s, ^s, and ^l- The first two are 
expressed as^^^'^^'^ 

1 „- V3 . 



cos 9s (t) 



■ cos t/s 



■ sin 6*3 cos[(/)(t) - 0s] 



TT _j / \/3cos6's + sin6lscos[(/)(t) - 0s] 

Psit) = ai H h tan = ^- 

' 12 I 2sine'ssin[0(t)-0s] 



(A.l) 
(A.2) 



For the polarization angle ipsit) (see Eq. (85)), first z ■ N = cos6's- Next when we 
neglect the spin precession effects, L is constant and L ■ L ■ N, and N ■ {L x z) 



are given as 



169,223 



(A.3) 
(A.4) 



i • 2; = i cos6'l - ^ sin6'L cos[0(t) - 0l], 
L ■ N ~ cos^L cos^s + sin^L sin^g cos(0l ~ 0s): 
N ■ {L X z) ^ isin^Lsin0ssin(0L - 0s) 
-\/3 

— — cos0(t)(cos0L sin 0s sin0g — cos 0s sin^L sin0L) 
■\/3 

— — sin 0(t) (cos 6*3 sin^L cos0l — cos^l sin 6*3 cos0g). (A. 5) 



Appendix A.2. Precessing Binaries 

The precession equations for circular orbit binaries are^^'^ 



1 



4mi + 3m2 4m2 + 3mi 

"Jl H n '->2 



2mi 



2m2 



3 1 
2a 



^,{S2-L)Si + iSi-L)S2] xL 



X L 

32/^ /Mt 
5 a \ a 



5/2 



Si 
S2 



1 
1 



4mi + 3to2 

2m 1 
4m2 + 3toi 

2m2 



(/iVMtai) 



X Si 



X S. 



1 



i52-^(52-i)X 
lSi-^-{Si-L)L 



(A.6) 
X 5i,(A.7) 
X 52.(A.8) 
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The first term of each equation represents the spin-orbit interactions and the second 
term represents the spin-spin interactions. The last term of Eq. (A. 6) is the angu- 
lar momentum loss due to the radiation reaction. This changes the total angular 
momentum J = L + Si + S2 as 

a \ a J 

Prom Eq. (Ill), the quantities {L ■ TV), {L ■ z) and [N ■ {L x 2;)], which are 
needed to compute the polarization angle -^g [t) in the beam-pattern coefficients 
and , are expressed as^^^ 

f - 17 '■\ \ , 1 - 2(Jo • i:)cos^j . 

L ■ z = [Jq ■ z) cos Al H = sm Al cos a 

2 sin 6'j 

+ ^ J sin Al sin a, (A.IO) 

sin Oj 

f jvr tT lin \ , cos^s - (Jo • -/V)cos0,j . 

L ■ N = (Jo • N) cos Al H = sm Al cos a 

sin 0} 

(Jo X -iV . ^ . 

H = smALsma, (A. 11) 

sin 0j 

/i- ~x /V X N ■ {% X z) - N ■ [Jo X z)cos9j . ^ 

N ■ [L X z) ^ N ■ (Jo X z) cos Al H = sm Al cos a 

sin 6] 

N ■ (Jo X z) X z , ^ , / . , „x 

+ ^ - ^ sin Al sin a, (A.12) 



sm Uj 



where 



1 - a/3 

Jo • 2 = - cos 6'j —smejcos[(j){t) - (A. 13) 

Jo ■ N = cos 6'j cos6's -I- sin^j sin6's cos((^j — (fs), (^-14) 
N ■ {Jo X z) ~ sin^j sin6's sin((^j — <j>s) 

— — cos (f>{t) (cos 9 J sin 6'g sin 0s — cos 9s sin 6'j sin 0j) 
•\/3 

— — sin 0(<) (cos 6'g sinflj cos0j — cos0j sin 6'g cos (/)s),(A.15) 
(Jo X ■ N — sin^g sin^j sin((^j — 0s), (A. 16) 

iV.(lx £) = ^sin^ssin(0(O-0g), (A.17) 
/3 - - 

(Jo X I) • 2 = ^ sin^j sin(0(i) - 0j), (A. 18) 

AT • (Jo X z) X z = — i sin^j[\/3cos^g cos{(f>{t) — 0j} -I- sinSg cos(0j — 0g)] . 

(A.19) 
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When spins are zero, sinAL = 0, = and <t>.^ ~ 0l- Then, Eqs. (A. 10), (A.ll) 
and (A. 12) each reduces to Eqs. (A. 3), (A. 4) and (A. 5), respectively. 
L ■ u and TV • (i x it) in Eq. (127) are given as^^^ 

L-u = -N ■ {L X Ja) (A.20) 

sinAL /" " NxSinA^ 

= AT • (Jo X zj— — — cosa — (cos^s — cos0j(Jo • N))— — — sina, 
sm Oj sm 

N-{Lxu)^ cosAl - (£ • iV) • (Jo • TV), (A.21) 
where TV • (Jo x z) and Jo • TV are given as Eqs. (A. 15) and (A. 14), respectively. 

Appendix A. 3. GW Energy- Momentum Tensor in BD Theory 

The energy-momentum tensor for GWs, t^^'^ , in BD theory is given as"'^^'^ 

+ ^XX''^'' + ^X'^X^'' + ^xDx'^^ + X.aX^'^'" - X,oX"^^'^^ 

-x^.^x'^"''' + x^,;3x'''''" - x.o'^^x-'^^" + ^x'^^Dx - 2xL"nx''^" 

+ \ri^'' (^Xa/3,^X"'''^ - \X,C.X^ + 2Xa/3nx"^ - Xa^,7X"^'^ " X^X + Xa^X"") 

+^vUx^- + U^x^'' - 2^.„'(^x"^" - 2(^,a(x"^^'''^ - X^"") + >,a0X"'' 

+ (2WBD - l)^'^ip-'' - ^W'""" + J?^" 

+0(x',xV,X^',V'') (A.22) 
with ip = 4>I4>Q- 
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